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Predictions for charged hadron, identified fight hadron, quarkonium, photon, jet and 
gauge bosons in p+Pb coUisions at y^s^Tjv" ~ ^ TeV are compiled and compared. When 
test run data are available, they are compared to the model predictions. 
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12.38.Bx, 25.75.Bh, 25.75.Cj, 13.87.-a 

1. Introduction 

Here predictions for the upcoming p+Pb LHC run at y/s^ — 5 TeV compiled 
by members and friends of the JET ConaboratiorP are presented. The test run 
data pubhshed by the ALICE CoUaboratiorP'^ are compared to model calculations 
available before the test run. Most calculations are for midrapidity and minimum 
bias collisions. Other results at different rapidities and centralities are presented 
when available. 

The predictions presented here, as well as the corresponding discussion, were 
made assuming that the proton circulated toward positive rapidity and the nucleus 
toward negative rapidity, p+Ph collisions, similar to fixed-target configurations. In 
the fixed-target configuration, the low x nuclear parton distributions are probed 
at positive rapidity. Throughout this paper, many of the results shown have been 
adjusted to the Ph+p convention of the ALICE data, as described in Ref.,^ with 
the nuclear parton density probed at high x at forward (positive) rapidity and low 
X at backward (negative) rapidity. The cases where the results still appear with 
the assumption that the proton circulates in the direction of positive rapidity are 
explicitly noted. 

This paper is organized in the following fashion. Section [2] describes the mod- 
els that specifically address charged particle production. These include saturation 
approaches, event generators, and perturbative QCD-based calculations. Section [3] 
compares results obtained from models described in Sec. [2] with each other and 
with the available data. The next several sections present predictions for specific 
observables including identified light hadrons (Sec.|4|), quarkonium (Sec. [5]), direct 
photons (Sec. [6]), jets (Sec. [7]), and gauge bosons (Sec. H]). 

2. Model descriptions 

In this section, the models used to obtain the results for charged particle distribu- 
tions, dNch/di], dNch/dpT and the nuclear suppression factor Rppb as a function 
of pt and ij. The first calculations described Sec. 12.11 are saturation or Color Glass 
Condensate (CGC) based. The next set of calculations are event-generator based 
with results from HIJING in Sec. [231 HIJINGBB2.0 in Sec. EH and AMPT in Sec. [231 



5 
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Finally, calculations based on coUinear factorization in perturbative QCD are de- 
scribed in Sees. 12.6] and 12.71 

In the following Section, Sec. |3l the predictions will be compared to the ALICE 
Ph+p test run data^ ^ in September 2012. 

2.1. Inclusive hadron production in the rcBK-CGC approach (A. 
Rezaeian) 

In the Color Class Condensate (CGC) approach, gluon jet production in p + A 
collisions can be described by fcT-factorization,'^ 

where Cp — {N^ — l)/2Nc, Nc is the number of colors, xi^2 — (pT/\/s)e^^, pr 
and y are the transverse momentum and rapidity of the produced gluon jet, and 
^/s is the nucleon-nucleon center-of-mass energy. The unintegrated gluon density, 
(j)'^{xi; kr) , denotes the probability to find a gluon that carries fractional energy 
Xi and transverse momentum kx in the projectile (or target) A. The unintegrated 
gluon density is related to the color dipole forward scattering amplitude, 

4>A kr) = / d^br d^VT e'^^'"'^ V^Ma (a:.; ^t; &t) , (2) 

V J as [2itY J 

with 

Ma [xi] rr; br) = 2A/'f (xi; rr] &t) - M3 i^t'^ &t) , (3) 

where tt is the transverse size of the dipole and &r is the impact parameter of the 
collision. The subscript T stands for the transverse component. The dipole scatter- 
ing amplitude Mf satisfies the nonlinear small- cc JIMWLK evolution equations ,1^^^^ 
see below. 

In the fcy-factorized approach, partons in both the projectile and target are as- 
sumed to be at very small x so that the CGC formalism and small-x resummation 
is applicable to both the projectile and the target, assuming the projectile pro- 
ton moves in the direction of forward rapidity. This approach is valid away from 
the projectile fragmentation region. However, to treat the projectile fragmentation 
region in the forward region, an alternative approach developed in Refs.,'^'^the so- 
called hybrid approach, is better suited. In this approach, the projectile is treated 
perturbatively within the standard coUinear factorization scheme using the stan- 
dard DGLAP picture while the target is treated employing CGC methods. The 
cross section for single inclusive hadron production at leading twist in asymmetric 
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collisions such as p + A in the CGC approach is given by 



drjd'^PT 



K 



1 dz 



Pt, 



Xifg{xi,flp)NA{x2, —)Dh/g{z, ^Fi-) 



Pt, 



+ T.qXlfq{xi,fIp)NFix2, —)Dh/g{z, ^p,-) 



27r2 



dz z 



2 4 

xp z Pt Jk. 



d^kTk'^NF{kT:X2) I 



d^ 



,Xi 



(4) 



A if-factor has been introduced to effectively incorporate higher-order corrections. 
The parton distribution function of the proton, fj(x,y?p), depends on the light- 
cone momentum fractions x and the hard factorization scale ^.p- The function 
D / ^{z , ^-pr) is the fragmentation function (FF) of parton i to become final-state 
hadron h carrying a fraction z of the parent parton momentum at fragmentation 
scale /^Fr- The inelastic weight functions, Wi/j, and the DGLAP splitting functions, 



i/j, are given in Refl^The longitudinal momentum fractions xi and X2 are 
Pt 



Xf 



Xl 



Xf 
z 



X2 — xie 



-2», 



(5) 



The strong coupling multiplying the inelastic term in Eq. (j4|), is denoted a™. 
The superscript "in" is employed to differentiate it from the running value of 
in the rcBK equation. In the hybrid formulation, the strong coupling in the dilute 
regime of the projectile can differ from that in the rcBK description of the dense 
target (or dipole scattering amplitude) . The scale at which a'/ should be evaluated 
cannot be determined in the current approximation, a full NNLO calculation is 
required. The effects of different choices of a^" will be considered later. 

In Eq. (U), the factorization scale is assumed to be the same in the frag- 
mentation functions (/xpr — fip) and the parton densities. In order to investigate 
the uncertainties associated with choice oi ^f, several values oi fiF are considered: 
Hf = '^Pt; Pt; and pt/2. 

In Eq. the amplitude Np (Na) is the two-dimensional Fourier transform 
of the imaginary part of the forward (proton direction) dipole-target scattering 
amplitude, AfA{F), in the fundamental (F) or adjoint (A) representations. 



NA{F){x,kT) ^ d fe 



-ikT - T 



1 



-A/-^(^)(r,r^ln(^)) 



(6) 



where r = |?^ is the dipole transverse size. The dipole scattering amplitude Ma{F) 
incorporates small-x dynamics and can be calculated using the JIMWLK evolution 
equation.l^'IS] jj]^ the large Nc limit, the coupled JIMWLK equations are simplified to 
the Balitsky-Kovchegov (BK) equation,EHlll g, closed- form equation for the rapidity 
evolution of the dipole amplitude. While a numerical solution of the full next-to- 
leading logarithmic expressions is not yet available, the running coupling corrections 
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to the leading log kernel, the so-called running-coupling BK (rcBK) equation has 
been very successful in phenomenological applications The rcBK equation has 
the following simple formEiHlllEll 

^dHxl'/xt " / ^™"(^^'^^i'^^2) [A/'A(F)(ri,x) +A/-^(^)(r2,x) 

- J^A(F) {r,x) - Ua(f) {ti , x) Ua(f) {r2 ,x)] , (7) 

where f2 = r — ri. The rcBK equation only describes the rapidity /energy evolution 
of the dipole, the initial profile and parameters of the dipole still need to be mod- 
eled and constrained by experimental data. The initial condition for the evolution 
generally takes a form motivated by the McLerran-Venugopalan model,'^^'^ 



N{r,Y^O) = 1-exp 



4 V Ar 



(8) 



where the onset of small- a; evolution is assumed to be at xq = 0.01, and the in- 
frared scale is A = 0.241 GeV.'^ The only free parameters are 7 and the initial 
saturation scale Qos, with s ~ p and A for proton and nuclear targets, respectively. 
Unfortunately, the current global set of small-a: data are very limited and thus can- 
not uniquely fix the initial dipole parameters.'^ This problem is more severe for 
determining the dipole scattering amplitude on nuclear targets, leading to rather 
large unavoidable theoretical uncertainties on CGC predictions ioT p + A collisions 
at the LHC. In Ref.jl^flJ a simple scheme to test the CGC dynamics at the LHC was 
proposed. This scheme will be used to calculate the results shown later on. 



2.2. IP-Sat (P. Tribedy and R. Venugopalan) 

The impact parameter dependent dipole saturation model (IP-Sat is a refinement 
of the Golec-Biernat-Wusthoff dipole modeP^'^to give the right perturbative limit 
when the dipole radius tt 0^ It is equivalent to the expression derived in the 
classical effective theory of the CGC, to leading logarithmic accuracy .'^S'^H The 
proton dipole cross section in this model is expressed as 



(Phr 



rT,x,bT) = 2 



1 — exp 



2Nc 



Tj^asipL )xg{x,fi )Tp{bi 



Here the scale is related to dipole radius tt as 



(9) 



(10) 



where the leading order expression for the running coupling is 

2, 127r 



(11) 



(33-2n/)log(MVA|cD) 

with Uf = Z and Aqcd = 0.2 GeV. The model includes saturation as eikonalized 
power corrections to the DGLAP leading-twist expression and may be valid in the 
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Table 1. Parameters of the IP-Sat model 
obtained from fits to HERA datal^H 

Bp (GeV-^) fiQ (GeV^) Ag Ag 
4.0 1.17 2.55 0.020 



regime where logs in dominate logs in x. For each value of the dipole radius, the 
gluon density xg{x, ^i^) is evolved from fi^ to ^j? using the LO DGLAP evolution 
equation without quarks, 

1 



d log /i^ 



27r 



(12) 



Here the gluon splitting function is 

z 1 — z 



1-z) 



+ z{l-z) 



The initial gluon density at the scale /ig is taken to be of the form 

xgix,fil)^AgX-^^{l-xf-^ . 



(13) 



(14) 



An important feature of the IP-Sat model is the 6-dependence of the dipole cross 
section, introduced through the gluon density profile function Tpibx)- This profile 
function, normalized to unity, is chosen to have the Gaussian form 



Tpibr) 



1 



■ exp 



f 



(15) 



2t:Bg V2Sg, 

where Bq is a parameter fit to the HERA diffractive data. It corresponds to (fe^) — 
2Bq, the average squared gluonic radius of the proton. 

The IP-Sat model parameters are obtained from optimal fits to HERA data.^^ 
The parameters used in this work are listed in Table 12.21 The parameters of the 
initial gluon distribution are determined from fits to the HERA F2 datcPSI122l with 

^ 1. The value of Bq is determined primarily from the J/-^ ^-distributions 
measured by ZEU^^SI and HI. 1^ With these parameters, excellent agreement with 
the HERA exclusive vector meson and DVCS data is obtained. For a detailed 
comparison of this model to the HERA data, see Ref.'22l A more recent fit to the 
combined ZEUS and HI data has been performed in Ref.^^^ 

The IP-Sat model successfully describes the bulk features of the p+p and A + A 
data over a wide range of center-of-mass energies from RHIC to LHC as well as the 
features of the RHIC d+Au data.'^'^ It also provides the basis for the IP-Glasma 
modeP^'^ of initial conditions in heavy-ion collisions. 

In the IP-Sat model, the dipole-nucleus cross section in a large nucleus can be 
approximated as 



d2 



ST 



1 — exp 



ATa{st) 



'dip 



{rT,x) 



(16) 
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where ATji{sT) is the transverse density of nucleons inside a nucleus and cr^ip(''Tj 
is obtained by integrating the dipole-proton cross section in Eq. (j9|) over the impact 
parameter distribution in the proton. This form of the dipole-nucleus cross section 
was previously shown to give reasonable fits to the limited available inclusive fixed- 
target e + A data.tSS 

111 p + A collisions, the LO inclusive gluon distribution can be expressed a^^ 

dNP'^ibr) _ Aas 1 f d^kr f ,2 d(l>p{xi,kT\sT) d(f>A{x2 

« ■5T yr yr -(l') 



dy d^PT ttCf Pt J (Sti")^ 7 d'^sx d^ST 

This equation is a generalization of the well known /c^-factorized expression for 
inclusive gluon productiorf^ to include the impact parameter dependence of the 
unintegrated gluon distributions. Here Cp = {N^ — i)/2Nc is the Casimir for the 
fundamental representation. Using a relation between quark and gluon dipole am- 
plitudes strictly valid in the large Nc limit, the unintegrated gluon distribution 
in protons and nuclei can be expressed in terms of the corresponding dipole cross 
section measured in DIS as**^ 

.,A -.2 



dr'^{x,kT\sT) _ J d^rre^'^T.rr . 1 



d'^ST 



2 d^ST 



(18) 



2.3. HIJING2.1 (R. Xu, W.-T. Deng and X.-N. Wang) 

The HIJINCP^l'^ Monte Carlo is based on a two-component model of hadron pro- 
duction in high-energy p + p, p + A, and A + A collisions. The soft and hard 
components are separated by a cutoff momentum po in the transverse momentum 
exchange. Hard parton scatterings with px > Po are assumed to be described by 
perturbative QCD (pQCD), while soft interactions are approximated by string exci- 
tations with an effective cross section CTgoft- Inp+A collisions, the single jet inclusive 
cross section is proportional to the nuclear parton densities fa/A{x2,PTj ^); 

^ Ja^ = AW dy2d'^bTA{b) ^ xifa/p{xi,p'^)x2fa/A{x2,PT^b) — . 

a,b,c 

Here, xi,2 = prie^^^ + e^^^)/\/s are the fractional momenta of the initial partons 
while are the rapidities of the final parton jets. Higher order corrections are 
absorbed into the K factor. The nuclear thickness function is normalized to A, 
JdHTAib) = A. 

There are several cold nuclear matter effects that are considered in HIJING. The 
first is the shadowing effect. HIJING2.Cp^'^ employs a factorized form of the parton 
densities in nuclei,'^^ 

fa/A{x,fJ.%,b) ^ Sa/A{x, fJ.%,b)fa/A{x, fip) (19) 

where Sa/A{x, y?p,b) is the impact-parameter dependent nuclear modification fac- 
tor. However, the shadowing employed in HIJING2.0 does not include any ii\ de- 
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pendence as in e.g. Ref.l^ Therefore, shadowing effects in HIJING should only be 
valid at low and disappear at larger . 

The second cold matter effect included is the Cronin effect the enhancement 
of intermediate pt hadron spectra hi p -\- A collisions. Multiple scattering inside a 
nucleus can lead to the transverse momentum {kr) broadening of both the initial- 
and final-state partons. A fcT—kick is imparted to both the initial and final-state 
hard scattered partons in each binary nucleon-nucleon scattering. The fcr-kick of 
each scattering follows a Gaussian distribution. Fits to the fixed-target p + A data 
lead to an energy dependence of the Gaussian width, 

(4) = [0.141og(\/i/GeV) - 0.43] GeVVc^ • (20) 

This fcx-kick influences the final-state hadron rapidity distribution. After tuning 
the gluon shadowing parameter sj^ in HIJING2.1, the charged particle rapidity 
distribution, dNch/dr], in d+Au collisions at ^s„„ = 200 GeV can be described. 
The prediction for LHC energies can be obtained by extrapolation.^*^ 

In the default HIJING setting, p + A and A + A collision are decomposed into 
independent and sequential nucleon-nucleon collisions. Within each nucleon-nucleon 
collision, hard collisions are simulated first, followed by soft collisions. However, 
since the time scale for hard scattering is much shorter than soft interactions, such 
a sequence of hard and soft interactions within each binary collision might not be 
physical. In a revised scheme denoted by DHC (decoherent hard scattering), all 
the hard interactions in a p + A event are simulated first. They are subsequently 
followed by the soft interactions. As a consequence, the energy available in each 
hard scattering is no longer restricted by soft interactions. 

An additional cold matter effect arises from valence quark number conservation 
in the proton. In p+A collisions, the projectile proton will suffer multiple scatterings 
within the target nucleus. For each binary nucleon-nucleon collision, there is a finite 
probability for independent hard parton scattering involving initial partons from 
the projectile and target nucleons. Flavor conservation limits the availability of 
valence quarks from the projectile for each of these hard interactions. This effect 
can change the relative flavor composition of produced partons per average binary 
nucleon-nucleon collision. Since the gluon fragmentation functions are softer than 
those of the quarks, the increased fraction of produced gluon jets mp + A collisions 
can lead to suppression of the final-state high-p^ hadron spectra. 

Finally, jet fragmentation can also modify the final hadron spectra in p + A 
collisions. In the default HIJING setup, jet shower partons from initial- and final- 
state radiation are ordered in rapidity. Gluons are connected to the valence quark 
and diquark of the projectile or target nucleons as kinks to form systems of strings. 
These strings fragment into final-state hadrons using the Lund string fragmentation 
model.'^ In p + A collisions, the projectile can undergo multiple scatterings since 
its string systems have many more gluons attached to them than in p + p collisions. 
Hadrons produced by the fragmentation of such string systems are softer than those 
resulting from independent fragmentation of individual gluons. 
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2.4. HIJINGBB2.0 (G. G. Barnafoldi, J. Barette, M. Gyulassy, P. 
Levai, M. Petrovici, and V. Topor Pop) 

Monte Carlo models such as HIJINGl.oj^^'Sl] HIJING2.cPI3Sl and HI JINGBB2.(PHI1] 

have been developed to study hadron production in p + p, p + A and A + A colh- 
sions. They are essentially two-component models which describe the production of 
hard parton jets and the soft interaction between nucleon remnants. Hard jet pro- 
duction is calculated employing collinearly-factorized multiple minijet production 
within pQCD. A transverse momentum cut-off, pq, on the final-state jet production 
is introduced so that for pT < po the interaction is nonperturbative and is charac- 
terized by a finite soft parton cross section (Tsoft- The jet cross sections depend on 
the parton distribution functions parameterized from global fits to data!^!^ 

Nucleon remnants interact via soft gluon exchanges described by the string mod- 
el^^^t^and constrained from lower energy + e~ , +p, and p+p data. The hard 
jet pairs and the two excited nucleon remnants are connected by independent strings 
which fragment to resonances that decay to the final-state hadrons. Longitudinal 
beam-jet string fragmentation depends strongly on the values of the string tensions 
that control the quark-antiquark (qq) and diquark-antidiquark {qqqq) creation rates 
and strangeness suppression factors (7s). 

In HIJINGl.O and HIJING2.0, a constant (vacuum value) for the effective value 
of string tension, kq = 1.0 GeV/fm, is used. At high initial energy density, the 
novel nuclear physics is due to the possibility of overlapping multiple longitudinal 
flux tubes leading to strong longitudinal color field (SCF) effects. These effects are 
modeled in HIJINGBB2.0 by varying the effective string tension. SCFs also mod- 
ify the fragmentation processes, resulting in an increase of (strange)baryons which 
play an important role in the description of the baryon to meson anomaly. In order 
to describe the p + p and central Pb-|-Pb data at the LHC, we have shown that 
the energy and mass dependence of the mean value of the string tension should be 
taken into account .'^'^ Moreover, to better describe the baryon to meson anomaly 
seen in the data, a specific implementation of JJ loops, has to be introduced. For 
a detailed discussion, see Refs.'^'^'^ Similar results can be obtained by includ- 
ing extra diquark-antidiquark production channels from the strong coherent fields 
formed in heavy-ion collisions.'^ 

All HIJING-type models implement nuclear effects such as modification of 
the parton distribution functions, shadowing^ and jet quenching via medium- 
induced parton splitting. (CoUisional energy loss is neglected.*'^ In HIJINGl.O 
and HIJINGBB2.0, the Duke-Owen (DO) parameterizations of the proton parton 
densitie^^ is used to calculate the jet production cross section with pT > po- 
In both codes, a constant cutoff, po = 2 GeV/c, and a soft parton cross sec- 
tion, (Tsoft — 54 mb, fit the experimental p + p data. However, for A + A colli- 
sions in HIJINGBB2.0, an energy and mass dependence of the cut-off parameter, 
Po{s,A) = 0.416y(s)° i9Mo i28 GeV/c, was introducecPl^lSl] at RHIC and LHC 
energies in order not to violate the geometrical limit for minijets production per 
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unit transverse area. The p -\- p cutoff was kept constant at po = 2 GeV/c. In 
HIJING2.o!^E6la subsequent version of HIJINGLCpMsl the GRV parameterization 
of the proton parton densitied^ is implemented. The GRV small x gluon density is 
much higher than that of the DO parameterization. Here also an energy-dependent 
cutoff pois) and soft cross section cr soft{s) are also assumed in order to better de- 
scribe the Pb-I-Pb data at the LHC. The cutoff used in HIJINGBB2.0 varies from 
po = 1.5 GeV/c at the CERN SPS, = 20 GeV, to 4.2 GeV/c at = 5.5 TeV 
while that in HIJING2.0, with a more complex energy dependence, varies from 1.7 
to 3.5 GeV/c in the same energy range. 

One of the main uncertainties in the calculation of the charged particle multi- 
plicity density in Pb-|-Pb collisions is the nuclear modification of parton distribution 
functions, especially gluon distributions at small x. In HIJING-type models, the par- 
ton distributions per nucleon in a nucleus, fa/Ai^, A*!-), are factorizable into parton 
distributions in a nucleon, fa/N, and the shadowing function for parton a, Sa/A, as 
in Eq. (fT9| . The shadowing parameterization in HIJINGl.Cp^'^ is employed, 

fa/A{x) 



Sa/A{x) 



A fa/N (x) 

1 + 1.191og^/^A [x^ - 1.2x^ + 0.21a; 

10.8 



1 



g-.Vo.oi ^ (21) 



log(A+l) 

assuming the same dependence for quarks and gluons. The evolution of 
Sa/Aix, fip) is neglected. The parameter, Sa, which determines the shadowing for 
X < 0.1, the region with the strongest nuclear dependence, is Sq = 0.1. For x > 0.1, 
the A dependence is rather weak. The parameterization in Eq. ((2T|) agrees with the 
X dependence of the quark structure function at small and medium cc.'^'^ Because 
the first part of Eq. ()21|) has a weak A dependence, impact parameter dependence 
is only included on part proportional to s^. The impact parameter dependence is 
given as 



5 / h 



where Ra is the radius of the nucleus and Sa = Sq = Sg = 0.1. 

The LHC Pb-|-Pb data at ./s^;^ = 2.76 TeV^^I^ indicate that impact-parameter 
dependent shadowing is required to understand the centrality dependence of the 
charged particle multiplicity density at midrapidity. These data place an indirect 
and model-dependent constraint on quark and gluon shadowing. Therefore, it is 
important to directly study quark and gluon shadowing in p + A collisions at the 
LHC. 

In contrast, in HIJING2.0,E3E£lthe factor (A^/^-l) is raised to the power 0.6 and 
a stronger impact-parameter dependence, different for quarks and gluons, Sq = 0.1 
and Sg — 0.22 — 0.23 respectively, is used to fit the LHC data. This stronger gluon 
shadowing requires jet quenching to be neglected.'^'^ 
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All HIJING-type models assume scale-independent shadowing (independent of 
Q^). This approximation could break down at sufficiently large scales due to the 
dominance of gluon emission in the DGLApl^ evolution equation. At Q = 2.0 and 
4.3 GeV/c, typical scales for minijet production at RHIC and LHC respectively, 
low X gluon shadowing varies by « 13% in the EPS09 LO parameterization.33 

2.5. AMPT (Z. Lin) 

The multiphase transport model AMPTfS3 was also used to calculate the yields and 
Pt spectra of particles produced in p + p and p+Pb collisions, as well as the nuclear 
modification factors i?ppb in p+Pb collisions. The flow coefficients have also been 
calculated. Both the default, AMPT — def and the string melting, AMPT — SM, versions 
of AMPlfS^ have been employed. 

In the default version of AMPT, AMPT — def, only minijet partons rescatter in the 
parton stage. After that, Lund string fragmentation is used for hadronization with 
the hadron cascade setting in at relatively high energy density. The cutoff time for 
the hadron cascade in these simulations is 30 fm/c (NTMAX = isd^. On the other 
hand, the string-melting version of AMPT, AMPT — SM, converts the usual initial-state 
hadronic strings to partonic matter when the energy density in the collision overlap 
volume is expected to be higher than that of the QCD phase transition. AMPT — SM 
also uses a simple quark coalescence model to describe bulk hadronization of the 
resultant partonic matter. Thus, secondary interactions are typically dominated 
by hadron interactions in AMPT — def while dominated by parton interactions in 
AMPT - SM. 

Using the default HIJING parameters for the Lund symmetric splitting function 
gives reasonable charged particle pseudorapidity distributions, dNch/drj at central 
values of pseudorapidity for Pb-l-Pb collisions at LHC energies.'SSl Therefore the 
same values (a = 0.5 and b = 0.9 GeV~^) are used for both p + p and p-|-Pb 
collisions, along with the same values of the strong coupling constant and parton 
cross section as in Ref.l^ 

In these calculations, p+p events are minimum-bias, including diffractive events. 
The MB p-|-Pb events include no restrictions on impact parameter. The nuclear 
modification factor, RpPh, as a function of px is obtained by dividing the p+Ph 
distribution by the p+p distribution, both calculated with the same version of 
AMPT, normalized by the number of binary N + N collisions, iVcoU- The number of 
collisions is assumed to be equal to the number of participant nucleons in the Pb 
nucleus (A^coii = A^part)- The collision centrality is defined according to the number 
of charged hadrons within \ri\ < 1. 

2.6. Leading-order pQCD calculations (Z.-B. Kang, I. Vitev, H. 
Xing) 

The details of the calculations described here can be found in Ref.^S^ A summary 
is given here. To leading order in the framework of factorized perturbative QCD, 
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single inclusive hadron production in p + p collisions, p{pi) + p{p2) ^ h{ph) + X, 
can be written a^^H 

da ^^oil v-^ f dxi ,9, „ , , 9 < /" dx2 ,9, „ , ,9 - 

J-^p / —'^''T^fa/Nixl,^J / — d^fcr^ /b/jv (x2 , fc^ ) 

X ^Dh/c{zc)Hab^cis,i,u)S{s + i+u), (23) 

where y and are the rapidity and transverse momentum of the produced hadron 
and J2abc ^^^^ parton flavors. In Eq. (|23)) . s = {pi + ^2)^; Dh/c{zc) is 
the fragmentation function (FF) of parton c into hadron ft,; Hab^c{s, t, u) are hard- 
scattering coefficient functions dependent on the partonic Mandelstam invariants 
s, £, u!^ A phenomenological K factor is included to account for higher-order QCD 
contributions. The parton distribution functions, ,fa.b/N{x,k'^), are dependent on 
the longitudinal momentum fraction x and the partonic transverse momentum ■ 
The fcT-dcpendence is included in order to incorporate the Cronin effect in p + ^ 
collisions. A Gaussian form is assumed, 1^ 

fa/N{xi,k^T,) = fa/N{xi)^^e-''-i/'''-\ (24) 

where fa/Ni^i) are the usual collinear PDFs in a nucleon. The factorization scale 
dependence has been suppressed in the arguments of fa/N- 

Inp + p collisions, (fc|.)pp = 1.8 GeY^/c^. The CTEQ6L1 PDFplare used with 
the fDSS parameterization of the parton-to-hadron fragmentation functions.!^ The 
factorization and renormalization scales are fixed to the transverse momentum of 
the produced particle, /ip = I^-r = Pt, and are suppressed in Eqs. and (1^^ . 
An 0(1) if -factor is found to give a good description of hadron production at both 
RHIC and LHC energies. 



2.6.1. Cold nuclear matter effects 

The p+A (e.g. d-|-Au or p+Ph) nuclear modification factor, RpA, is typically defined 
as: 

n -1 



RpA 



da, 



pA 



dyd^pj 



da, 



pp 



(25) 



_{Ncoii) dyd'^PT_ 

where (Ncow) is the average number of binary nucleon-nucleon collisions. The devi- 
ation of RpA from unity reveals the presence of cold nuclear matter (CNM) effects 
in p + A collisions. 

A variety of CNM effects can affect particle production. This section describes 
those that arise from the elastic, inelastic and coherent scattering of partons in 
large nuclei.'^ The proton and neutron composition of the interacting nuclei are 
also accounted for. In particular, these effects include isospin, the Cronin effect, 
cold nuclear matter energy loss and dynamical shadowing. These effects have been 
well documented in the literature. Their implementation is briefly described here. 
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Isospin The isospin effect can be easily accounted for on average in the nuclear 
PDFs for a nucleus with mass number A and charge Z by!^ 

fa/A{x) = ^fa/p{x) + " fa/n{x) , (26) 

assuming no modifications of the parton densities. In Eq. (|26|). fa/p{x) and fa/n{x) 
are the PDFs in a proton and a neutron, respectively. The neutron PDFs are related 
to those in the proton by isospin symmetry. 

Cronin effect The Cronin effect has been well documented.'^ It can be modeled 
by initial-state multiple parton scatterings in cold nuclei and the corresponding 
induced parton transverse momentum broadening.'^EIl particular, if the PDFs, 
fb I A {^2 , ) ' tiS'^s a normalized Gaussian form, the random elastic scattering in- 
duces further fc^-broadening in the nucleus: 

(fcT.)pA = (4.)pp + (|^)c. (27) 

Here A:t2 is the transverse momentum component of the parton prior to the hard 
scattering, C = ln(l + Sp"^), S = 0.14 (GeV/c)-^, fi^ = 0.12 (GeV/c)^, and Xg = 
{Cf/Ca)\ = 1 fm. These parameters describe the RHIC data reasonably well. 

Cold nuclear matter energy loss As the parton from the proton undergoes 
multiple scattering in the nucleus before the hard collision, it can lose energy due 
to medium-induced gluon bremsstrahlung. This effect can be easily implemented 
as a shift in the momentum fraction in the PDFs 

fq/pixi)-^ fq/p\Y^ ) ' fa/pi"^^) f9/p\T^ ) ■ 

Ideally, Eq. ([25)) should include a convolution over the probability of cold nuclear 
matter energy loss.^^S However, concurrent implementation of this distribution to- 
gether with the Cronin effect and coherent power corrections is computationally 
very demanding. The main effect of fluctuations due to multiple gluon emission 
is an effective reduction in the fractional energy loss Cq.gcff relative to the mean 
value (e^.g) = (^^^{AEi/E)) where the sum runs over all medium-induced gluons. 
Here ey.gcs = 0.7(e^.g). The average cold nuclear matter energy loss is obtained by 
integrating the initial-state medium-induced bremsstrahlung spectrum first derived 
in Ref.'^ It also depends on the typical transverse momentum transfer squared per 
interaction between the parton and the medium and the gluon mean- free path Xg. 
Therefore, the parameters are constrained to be the same as in the implementation 
of the Cronin effect, /i^ = 0.12 (GeV/c)^ and Xg — 1 fm. This calculation of initial- 
state cold nuclear matter energy loss has been shown to give a good description of 
the nuclear modification of Drell-Yan production in fixed-target experiments.'^ 

Dynamical shadowing Power-suppressed resummed coherent final-state scatter- 
ing of the struck partons leads to shadowing effects (suppression of the cross section 



January 23, 2013 2:55 WSPC/INSTRUCTION FILE predictionsall 



Predictions for p+Pb Collisions at = 5 TeV 17 



in the small- a; region)?^ The effect can be interpreted as the dynamical generation 
of parton mass in the background gluon field of the nucleus Thus 

x^x(^ + c /^^^'^.'^^ ^ , (29) 

where x is the parton momentum fraction in the lead ion, Cd — Cf{Ca) if the parton 
d = q{g) in the 2-^2 parton scattering ab cd, and represents the characteristic 
scale of the multiple scattering per nucleon. At RHIC energies, y/s^ = 200 GeV, 
— Cf/Ca^I — 0.12 GeV^f^ gives a good description of the nuclear modification 
in d+Au collisions for both single hadron and dihadron productionP^l 



2.7. Initial-state Shadowing ( G. G. Barnafoldi, J. Barette, M. 
Gyulassy, P. Levai, G. Papp and V. Topor Pop) 

The calculations in this section use the kTpQCD_v2.0 code, based on a 
phenomenologically-enhanced, perturbative QCD improved parton model described 
in detail in RefsP^"^ The main feature of this model is the phenomenologically- 
generalized parton distribution function employed to handle nonperturbative effects 
at relatively low-x and small px- The model includes intrinsic kx broadening with 
the average kx left as a free parameter to correct for nonperturbative effects. The 
kx value is determined from p + p data over a wide range of energies. Within the 
framework of this model, the fcr-broadening mp + A and A + A collisions is related 
to nuclear multiple scattering and can generate the Cronin enhancement!^'^ that 
appears within 3 < px < 9 GeV/c from SPS to RHIC energies. 



2.7.1. Theoretical Background 

The kTpClCD_v2.0 code calculates the invariant cross section for hadron production 
in p + p, p + A and A + A collisions at LO or NLO in the fcT-enhanced pQCD- 
improved parton model assuming collinear factorization. The code provides a Monte 
Carlo-based integration of the convolution,!^ written here for p + p collisions. 



d^px 



1 



■E 

abc 



-{l~V)/z^ 
VW/z^ ^(1 



dv 



— f dZr 



d'kx, 
da 



d^kx^ fa/p{xi 



VW/vz^ W 

kn, P%) fb/p{x2,kx2, f^p) 



(30) 



dv 



S{1 



W) ^ KabA^,V,W,pF, PR: PFi) 



Here da /dv represents the Born cross section of the partonic subprocess ah — >■ cd 
while Kab. c{s, V, w, ^iF, IJ-R, IJ-Fi) is the next order correction term. The proton and 
parton level NLO kinematic variables are (s, V, and W) and (s, v, and w) respec- 
tively.^^ 84 86 rpj^g various scales are /ip, the factorization scale; /i^, the renormal- 
ization scale; and /ipr the fragmentation scale. The factorization and the renormal- 
ization scales are related to the momentum of the intermediate jet, pf — PR =^ i^Pq 
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where k = 2/3, Pq — px/zc and Zc is the fraction of parton c momenta trans- 
ferred to the final hadron h. The fragmentation scale is related to the final hadron 
momentum by /ipr = i^Pt- 

The x-dependent proton parton distribution functions, fa/pix, iSp), defined in 
the infinite momentum frame, are generalized to three dimensions by incorporating 
an initial fcr dependence, 

fa/p{x, kr, fip) = g{kT)fa/A{x, Hp) ■ (31) 

The two-dimensional initial transverse momentum distribution, g^kx), with intrin- 
sic parton kr employed in these calculations is described in Refs jllllllll 87l[88| r^j^g 
distribution is described by a Gaussian, 

Here (fc|i) is the width of the kx distribution, related to the magnitude of the average 
parton transverse momentum by {k'^) = 4(fcT)^/7r. This treatment was successfully 
applied at LO in Ref . along with a iiTjct-based NLO calculation!^!^ In order to 
reproduce results for NN collisions at relatively low x, {k'^) = 2.5 GeV^/c^ was 
required. 

The LO or NLO fragmentation functions, (zc, Mpr); the probability for 
parton c to fragment into hadron h with momentum fraction Zc at fragmentation 
scale fiFi-. The MRST(cgj^ parton densities are used in Eq. (|3T|) . along with the 
KKP parameterization'*^ of the fragmentation functions. Both these sets can be 
applied at relatively small scales, = yitp^ « 1.25 GeV^. Thus the results obtained 
in these calculations are applicable for pr > 2 GeV/c. 

2.7.2. Incorporating Initial-State Nuclear Effects in p + A and A + A 
Collisions 

Proton-nucleus and nucleus-nucleus collisions can be described by incorporating 
the appropriate collision geometry and and nuclear shadowing. In the Glauber 
framework, the cross section for hadron production in an j4 -|- A' collision can be 
written as an integral over impact parameter b: 

eJ§^ = / ^^'^'^ ^-(1^^- ^D^'^ ^<-((fc^)M,(fe^)MO . (33) 

Here the nuclear thickness function, Ta(6) — J dz pA(b, z), employing the Woods- 
Saxon density distribution, is normalized so that J d^bTji(b) = A. 

The p + p cross section from Eq. pO[) includes increased fcx widths relative 
to p + p collisions, Eq. p2l) . as a consequence of multiple scattering in nuclei, see 
Eq. ([M)) . The increased width of the fcr distribution is taken into account by adding 
a function, hpA{b), that describes the number of effective N + N collisions at impact 
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parameter b, weighted by the average transverse momentum squared imparted by 
each collision, to the width in p + p collisions, {k'j.] 



pp, 



{ki)pA = {k^)pp + ChpAib) . (34) 

The function hpA{b) is expressed in terms of the number of collisions suffered by 
the incoming proton in the target nucleus, VA{b) = aNNTA^b), where aj^N is the 
inelastic N + N cross section: 

= ■ (35) 

[ zy„j — 1 otherwise 

For heavy nuclei, the maximum number of collisions is 3 < Vm ^ 4 with C — 
0.4 GeVVc2. 

Finally, the nuclear PDFs are modified by shadowing l^'^lHlll This effect, as 
well as the isospin asymmetry, are taken into account on average using the scale 
independent parameterization of Sa/Ai^) adopted from Ref.,'^ 



fa/A{x,tJ-%) = Sa/A{x) 



(36) 



where the neutron parton density, f a / n{x , IJ^\) , is related to that of the proton. 
Results are shown with the EKS98,^^ EPSOS^^ and HKN^- parameterizations, as 
well as with the updated HIJING parameterization:"^^ The EKS98, EPS08 and HKN 
parameterizations differ for quarks, antiquarks and gluons but are independent of 
impact parameter. The new HIJING parameterization differentiates between quarks 
and gluons but can include impact parameter dependence, as in Eq. (f22|) . 



3. Charged particles 

In this section, results on the charged particle multiplicity and pT distributions 
and the suppression factor Rpph as a function of pt are compiled. These results 
are compared with the ALICE test beam data where available. Other, related, 
predictions for charged particle observables are also shown. The upcoming p+Pb 
run at the LHC can place important constraints on models of the initial state. 

Note that the LHC magnet design requires the magnetic rigidity of the beams in 
the two rings to be the same. The proton beam, at 4 TeV, circulated in negative z- 
direction (toward negative rapidity) in the ALICE laboratory frame while a beam 
of fully-stripped Pb ions of (82/208) x 4 TeV/nucleon circulated in the positive 
z-direction (toward positive rapidity), implying Pb+p collisions if the first-named 
collision partner travels in the direction of positive rapidity, rather than p-|-Pb 
collisions, as the collisions are referred to throughout this text. This configuration 
resulted in a center of mass energy of ^s„„ = 5.02 TeV, moving with a rapidity 
difference, ^Vj^,^ = 0.465, in the direction of the proton beam. 

It is important to note that the most of the predictions shown here were orginally 
made assuming that the proton circulated toward positive rapidity and the nucleus 
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toward negative rapidity, similar to the fixed-target configuration and also the con- 
vention for d-f Au collisions at RHIC. In such cases, the low x nuclear parton dis- 
tributions are probed at positive rapidity. Here and throughout the remainder of 
this paper, the results shown have been adjusted to the Ph+p convention of the 
ALICE dateP unless otherwise explicitly noted. 

ALICE reported the primary charged particle pseudorapidity density in the 
laboratory frame, dA^ch/rf^yiab in non single-diffractive (NSD) p+Ph collisions. The 
lab frame pseudorapidity is defined as ryiab = — lntan(0/2) where 9 is the polar 
angle between the direction of the produced charged particle and the beam axis. 
The primary particles are due to both prompt production in the collision and strong 
decays. 

Calculations are typically performed in the center of mass frame but the pseu- 
dorapidity densities shown here have been also calculated in the lab frame. In the 
ALICE paper,'^ calculations in the center of mass frame were shifted by Ay„„ in 
the lab frame. This is only approximately correct since, at low px, the rapidity 
and pseudorapidity are not identical. The uncertainty on dNch/drjii^h due to the 
choice of frame was estimated to be less than 6%.^ It is worth noting that there is 
no ambiguity due to the calculational frame for identified particles thus the frame 
dependence is not discussed in later sections. 



3.1. Multiplicity distribution 

The calculations of the charged particle multiplicity distributions are described 
here. Saturation model predictions are discussed first, followed by event generator 
predictions and pQCD calculations with cold matter effects. 



3.1.1. Saturation Approaches (J. Albacete, A. Dumitru, H. Fujii, Y. Nara, 
A. Rezaeian, and R. Vogt) 

In the rcBK approach, used by Albacete and collaboratord^H ^nd Rezaeian for 
^pPbj'^'^the initial condition for the evolution of the dipole scattering amplitude 
can be written as (see also Eq. (jS])), 



M{r, a;o) = 1 — cxp 



1 ( 2 

r 

4 ^ 



^M^log(e+-l 



(37) 



There are three sets of unintegrated gluon distributions that are solutions of the 
rcBK small- a; evolution equations with the AAMQS initial condition in Eq. (|37|l . 
The values of 7 and Q^soi^o = 0.01) for protons that provide good fits to the e~ +p 
data are given in Table 13.1.11 Albacete et al. use all three initial conditions.'^ 
Rezaeian uses gl.llJpSll23 in the calculations of i?ppb in Sec. 13.31 but employs the 
b-CGC approaclPslllZl to calculate dNch/drj here. 



rcBK (J. Albacete, A. Dumitru, H. Fujii and Y. Nara) The calculation of 
dNch/dr] by Albacete et a/. employing 7 = 1.119, in both the center of mass and 



January 23, 2013 2:55 WSPC/INSTRUCTION FILE predictionsall 



Predictions for p+Pb Collisions at y/s^^ = 5 TeV 21 



Table 2. The AAMQS initial conditions used 
in the dipole evolution of the rcBK approach 
for protons 



Set 


Qloixo = 0.01) (GeVVc^) 


7 


MV 


0.200 


1 


gl.lOl 


0.157 


1.101 


gl-119 


0.168 


1.119 



lab frames, is shown in the dashed cyan curves in Fig. [TJ 




Fig. 1. Charged particle pseudorapidity distributions at ^s^^^ = 5.02 TeV in the CM (left) and 
lab (right) frames. The rcBK (dashed cyan) result is from Ref.lSH The band labeled Rezaeian 
(dot-dot-dash-dashed green with vertical bars outlining the uncertainty), described in Sec. 12.11 
is only shown in the center-of-mass frame. The band showing the IP-Sat result described in 
Sec. 12.21 is outlined by the dot-dot-dash-dashed black curves. The HIJING2.1 result without (NS, 
dot-dash-dash-dashed red) and with shadowing (sg = 0.28, solid red) and the HIJINGBB2.0 result 
without (dot dashed magenta) and with shadowing (dotted magenta) are also shown. Finally, the 
AMPT — def (dot-dash-dash-dashed blue) and AMPT — SM (dot-dot-dot-dash-dash-dashed blue) are 
given. The ALICE results from Ref.-^ are given on the right-hand side. The systematic uncertainties 
are shown in blue, the statistical uncertainties are too small to be visible on the scale of the plot. 



b-CGC (A. Rezaeian) The results of calculations by RezaeiarpSllMI in the center 
of mass frame are shown by the green curves in Fig. [TJ The vertical lines indicate 
the uncertainty in the predictions. The dipole forward (proton direction) ampli- 
tude is calculated in the b-CGC saturation model.-''' which has an explicit impact 
parameter dependence as well as approximately incorporates all known features of 
small-x physics. It describes the small-x HERA data, including diffractive data,!^ 
and also the RHIC and LHC data at small- cc.'^'^^ltifi^l xhis framework also provides 
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an excellent description of the charged hadron multiplicity in d+Au collisions at 
RHIC in addition to p + p and A + A collisions over a range of energiesP'® For other 
saturation model predictions in minimum bias p + A collisions at ^s„„ = 4.4 TeV, 
see Ref-EsH 

IP-Sat (P. Tribedy and R. Venugopalan) The calculation of the minimum- 
bias charged particle pseudorapidity distribution was performed in the kx- 
factorization approach using the IP-Sat modeP^ in Eq. (fT7|) . For the calculation 
of dNch/dr] in the lab frame, a constant rapidity shift of 0.46 has been applied. The 
bands shown in Fig. [T] arises due to uncertainties in the parameters. 

The charged particle multiplicity distributions obtained by the event generators 
HIJING, HIJINGBB2.0 and AMPT are now discussed. 

3.1.2. HIJING2.1 (R. Xu, W.-T. Deng and X.-N. Wang) 

The HIJING2.1 predictions in Fig. [T] are indicated by the dot-dash-dash-dashed 
curves without shadowing and the solid curves with Sg = 0.28. The difference 
between results with and without shadowing are largest for this calculation. 

3.1.3. HIJINGBB2.0 (G. G. Barnafoldi, J. Burette, M. Gyulassy, P. Levai, 
M.Petrovici and V. Topor Pop) 

The HIJINGBB2.0 predictions of dNch/di] in minimum bias collisions employ the 
values for the minijet cutoff and string-tension parameters of pq = 3.1 GeV/c 
and K = 2.0 GeV/fm. These values are determined from fits to p + p and A + A 
systematics from RHIC to the LHC, see 

R_gfP|54lf5l for details. Note that these 

calculations assume no jet quenching. 

The absolute normalization of dNch/drj is sensitive to the low px, Pt < GeV/c, 
nonperturbative hadronization dynamics encapsulated in the LuncP^'SH jETSElf^ 
string fragmentation constrained by lower energy + e~ , + p, and p + p data. 
The default HIJINGl.O shadowing parametrization leads to substantial reduction of 
the global multiplicity at the LHC. The HIJINGBB2.0 results without shadowing are 
substantially reduced relative to the same predictions with HIJINGl.cf^'^ because 
both the default minijet cutoff po ^ 2 GeV/c and vacuum string tension kq — 
1 GeV/fm used in HIJINGl.O are generalized to vary monotonically with y/s^ and 
A. As discussed in Ref. j'^'^'^ systematics of multi-particle production m A + A 
collisions from RHIC to the LHC are used to fix the y/s^ and A dependence 
of po and k. The resulting dependencies, po{s,A) = OAW^/^'^^^A'^-^'^^ GeV/c 
and k{s,A) = Ko(s/sO)° °4^'' i67 GeV/fnPl lead to po = 3.1 GeV/c and k = 
2.1 GeV/fm in p+Ph colhsions at 5.02 TeV. 

Constant values of the cutoff are employed in p+p collisions, independent of the 
incident energy, pg^ = 2 GeV/c, and string tension, Kpp — 1.9 GeV/fm. Note that, 
even without shadowing, the increase of po from 2 GeV/c in p + p collisions to 3.1 
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GeV/c in p+Pb collisions causes a reduction in the minijet cross section and hence 
the final pion (charged particle) multiplicity. Such a reduction is also required to 
fit the slow growth (by a factor of 2.2) in the A + A charged particle multiplicity 
from RHIC to LHCli^ll This reduction could be interpreted as phenomenological 
evidence for gluon saturation beyond leading twist shadowing. 

3.1.4. AMPT (Z.-W. Lin) 

The AMPT default and string melting calculations are shown in the blue curves in 
Fig. [T] The differences in the two scenarios is not large. Indeed it is much smaller 
than models with and without shadowing. These differences can arise from several 
sources: the relative rescattering strengths in the parton and hadron stages; the 
hadronization models; or a combination of rescattering and hadronization. 

3.1.5. Forward/Backward difference (R. Vogt) 

The event generator calculations produce distributions that do not show a strong 
forward/backward difference between the lead and proton peaks. The CGC-based 
calculations, however, show a much stronger dependence of the results on pseudo- 
rapidity, in both frames. 

The charged particle multiplicity results can be further quantified by comparing 
the measured to predicted particle density at midrapidity, near the proton peak, 
f?iab — —2, and the lead peak, r/iab = 2. The absolute values of dNch/drjii^h at 
Viah = ~2, and 2, along with the ratio R of the multiplicities at T^iab = 2 to 
??iab = —2) are given in Table [XT751 The ALICE results as well as the model results 
included both here and in the ALICE papei^^l are given. There are two tabulated 
values for both HIJINGBB2.0 NS (no shadowing) and WS (with shadowing). The first 
were shifted from the center-of-mass to the lab frame by the ALICE collaboration 
and the second were direct lab frame calculations. 

3.1.6. Centrality Dependence ofdNch/drj 

Here we present two calculations of the centrality dependence of the charged particle 
multiplicity distributions. 

AMPT (Z. Lin) Figure [2] shows the pseudorapidity rapidity distributions calculated 
in the center-of-mass frame by AMPT — def . Results are shown for p+p and minimum 
bias p-|-Pb collisions as well as for six different centrality classes. In AMPT, the p + p 
events are minimum bias events, including diffractive events. 

Since the ratio Rpph is normalized by a factor proportional to the number of bi- 
nary nucleon-nucleon collisions, A'coii, the number of collisions in the same centrality 
classes as included in Fig. [2] is shown in Table 13.1.61 In AMPT, A'coU is assumed to 
be equivalent to the total number of participant nucleons in the Pb nucleus, Np^^.^.. 
Because there is some model dependence in the definition of centrality bins, it is 
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Table 3. Comparison of values of rfAf^h /i^»?iab S't Vla.h = ~2, 0, 2 and the ratio 
'i^ch/'^'?lab|i7iab=2/'^^ch/'^^labl»)iab=-2' denoted by R below. The tabulated IP-Sat re- 
sult is the average of the upper and lower limits depicted on the right-hand side of 
Fig. [T] The * on HI JINGBB2.0 indicates that the calculations have been shifted to the 
lab frame by the ALICE Collaboration while the f are results provided by V. Topor 
Pop et al.. Adapted from Ref.111 







dAfch/rfrjiab 




R 




-2 





2 




ALICE 


16.65 ±0.65 


17.24 ± 0.66 


19.81 ± 0.78 


1.19 ±0.05 


Saturation Models 










IP-Sat 


17.55 


20.55 


23.11 


1.32 


KLN 


15.96 


17.51 


22.02 


1.38 


rcBK 


14.27 


16.94 


22.51 


1.58 


HIJING-basod 










HIJING2.1 NS (no shad) 


23.58 


22.67 


24.96 


1.06 


HIJING2.1 WS {sg = 0.28) 


18.30 


17.49 


20.21 


1.10 


HIJINGBB2.0 NS* 


20.03 


19.68 


23.24 


1.16 


HIJINGBB2.0 NS^ 


16.84 


16.39 


19.68 


1.16 


HIJINGBB2.0 WS* 


12.97 


12.09 


15.16 


1.17 


HIJINGBB2.0 WSt 


13.98 


13.71 


16.73 


1.20 


AMPT 










AMPT - Def 


19.07 


18.56 


21.65 


1.14 


AMPT - SM 


18.14 


18.10 


20.84 


1.15 


DPMJET 


17.50 


17.61 


20.67 


1.18 



Table 4. Centrality classes of p+Pb events from AMPT — SM The centrality is determined 
from the number of charged hadrons within |r;| < 1 in the center-of-mass frame. 



Bin 


(b> (fm) 


fcmin (fm) 


&max (fm) 


-"part 


part, m 


(A^chd'?! < 1)) 


MB 


5.84 


0.0 


13.2 


7.51 


5.37 


36.1 


0-5% 


3.51 


0.0 


8.8 


15.70 


12.12 


99.9 


5-10% 


3.76 


0.0 


9.4 


14.19 


10.71 


79.6 


10-20% 


4.00 


0.0 


9.9 


12.93 


9.59 


66.5 


20-40% 


4.56 


0.0 


12.0 


10.70 


7.69 


49.6 


40-60% 


5.65 


0.1 


13.2 


7.30 


5.03 


31.4 


60-80% 


7.08 


0.1 


13.2 


3.77 


2.50 


15.9 


80-100% 


8.08 


0.2 


13.2 


1.85 


1.11 


5.5 



worthwhile noting that other calculations of this same quantity may give somewhat 
different results. Here, the centrality of p-|-Pb collisions is defined according to the 
number of charged hadrons within \ri\ < 1. Table [37l.6l shows various information for 
each centrality class in the center-of-mass frame, including the average, minimum 
and maximum impact parameters, A^^a'rt: the number of participant Pb nucleons 
involved in inelastic collisions, N^^^.^ and the average number of charged parti- 
cles within |r/| < 1 from AMPT - SM. Note that the AMPT - SM values in Table l3T6l 
are essentially the same as those for AMPT — def except for the average number of 
charged particles (last column). 

In Fig. [21 the 40 — 60% centrality bin gives a distribution that is very close to 
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-10.0 -7.5 -5.0 -B.5 0.0 3.5 5.0 7.5 10.0 



Fig. 2. Charged particle pseudorapidity distributions in the center of mass frame for p + p, p+Ph 
minimum bias and six of the centraUty bins defined in Table [3^ 1.61 as calculated in AMPT. 

the min-bias result, in particular on the proton side. The most central bin is the 
most asymmetric distribution, as well as the largest in magnitude. On the other 
hand, the distribution for the 80— 100% centrality bin is symmetric around rjcm = 
and lower in magnitude than the p + p distribution. 

b-CGC (A. Rezaeian) In Fig. [3l the charged hadron multiplicity obtained from 
kr factorization in the b-CGC approach is shown the - 20%, 20 - 40%, 40 - 60%, 
60 — 80% centrality bins as well as minimum-bias collisions.'^ The impact parameter 
dependence of the saturation model is crucial for defining the collision centrality. 
The largest asymmetry in the multiplicity distribution is observed in more central 
collisions while, for peripheral collisions such as the 60 — 80% most central, the 
system becomes more similar to that produced in p + p collisions. This is reflected 
in the total charged hadron multiplicity distribution. 

Note that for all results shown in Fig. |3l a fixed minijet mass equal to current- 
quark mass is assumed for all energies/rapidities and centralities. Since the minijet 
mass is related to pre-hadronization/hadronization stage and cannot be obtained 
from saturation physics, it was fixed by fitting lower energy minimum-bias data. 
In very peripheral collisions where the system becomes more similar to symmetric 
p + p collisions, this assumption is less reliable. More importantly, one should also 
note that the kx factorization employed here is only proven in asymmetric p + A 
collisions at small x. Therefore, for more peripheral collisions, the current CGC 
prescription may be less reliable. 
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b-CGC [p+Pb, 5.02 TeV] 

50 
40 
30 
20 



10, 



-3 -2-10 12 3 

Fig. 3. Charged particle pseudorapidity distributions in the center of mass frame of p+Pb collisions 
at various centralities within the b-CGC saturation model. The ~ 5% theoretical uncertainties 
arising from fixing the overall normalization at RHIC are also shown. Prom 

3.2. Transverse Momentum distribution 

3.2.1. Compilation of midrapidity results (J. Albacete, A. Dumitru, H. Fujii, 
Y. Nara, R. Xu, W.-T. Deng, X.-N. Wang, G. G. Barnafdldi, J. 
Burette, M. Gyulassy, P. Levai, M. Petrovici, V. Topor Pop, Z. Lin 
and R. Vogt) 

The midrapidity, < 0.8, charged hadron pT distributions are shown on the 
left-hand side of Fig. H for rcBK,!^ HIJINGBB2.0, and AMPT. The HIJINGBB2.0 dis- 
tributions are similar to the rcBK results, albeit somewhat higher for px > 10 
GeV/c. The AMPT distributions, on the other hand, drop faster at low px than the 
other results but then become harder at high px- The AMPT results are essentially 
independent of whether string melting is included or not while the HIJINGBB2.0 
results without shadowing lie above those with shadowing. 

The right-hand side of Fig. |4] shows several options for cold matter effects 
in HIJING2.1 relative to p + p collisions. The p + p result is unsealed while the 
p-|-Pb curves with decoherent hard scatterings (DHC) without shadowing, DHC 
with shadowing, and shadowing only are separated from each other, starting from 
the p + p result, by a factor of 100. 

3.2.2. reBK at y = 0, 2 (J. Albacete, A. Dumitru, H. Fujii and Y. Nara) 

Fig. [5] shows the px spectrum in p + p (left) and p+Pb (right) collisions at different 
rapidities for the AAMQS initial conditions, see Table [3TTTJ Near central rapidity. 



1 1 1 1 1 1 1 1 1 1 
Min-bias 
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>« 20-40% tT^^ 
"40-60% T^0^ 
a-B 60-80% 
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rcBK band 7j = 
HIJING BB NS |7/|<0.8 ^ 
HIJING BB WS |7)|<0.B . 




Pt (GeV) (GeV) 



Fig. 4. (Left) Charged particle pT distributions at y/s^^ = 5.02 TeV. The sohd and dashed cyan 
curves outhne the rcBK band calculated by Albacete et aZ.I^Il The magenta curves, calculated 
with HIJINGBB2.0 as described in Sec. 12.41 are presented without (dot-dashed) and with (dotted) 
shadowing. The AMPT results, given by the dot-dash-dash-dashed (default) and dot-dot-dot-dashed 
(SM) blue curves, are described in Sec. 12.51 (Right) The charged hadron pj- distribution in p-l-Pb 
collisions with different HIJING2.1 options, scaled by the indicated factors to separate the curves. 
The p + p distribution is shown for comparison. 



fcji-factorization is employed while at forward rapidities (in the proton direction) 
the hybrid formalism is applied. The bands correspond to uncertainty estimates due 
to small variations in the scale entering the coupling and fragmentation functions. 




Fig. 5. The predicted transverse momentum spectrum in p -|- p (left) and minimum-bias p-|-Pb 
(right) collisions at y^s = 5 TeV at different rapidities (with the convention that the proton beam 
moves toward forward rapidity). From Ref.l^ 
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3.2.3. Flow Coefficients obtained with AMPT (Z. Lin) 

The flow coefficients w,i{2} {n — 2, 3, 4) in this study are calculated employing the 
two-particle cumulant method with v„ = ^ (cos — 0^)]), where cos \n{^(\)i — (f)j)] 
is averaged over all particle pairs in the specified phase space with both particles 
of each pair coming from the same event. 




p+p AMPT-SM 
p+p AMPT-det 
p+Pb MB AMPT-SM 
p+Pb MB AMPT-def 

_l , I , I , L 




p+p AMPT-SM 
p+p AMPT-def 
p+Pb MB AMPT-SM 
p+Pb MB AMPT-def 

J . \ . \ . L 




p+p AMPT-SM 
p+p AMPT-def 
p+Pb MB AMPT-SM 
p+Pb MB AMPT-def 



_L 



_L 



_L 





n 



Fig. 6. The flow coefBcients V2{2} (top left), V3{2} (top right), and V4{2} (bottom) of charged 
particles as a function of rj in p + p and minimum-bias p+Pb collisions. The results are shown for 
both AMPT - def and AMPT - SM. 



Figure [6] shows V2{2}, W3{2} and W4{2} for unidentified charged particles as a 
function of rj in p + p and minimum-bias p-|-Pb collisions from AMPT — def and 
AMPT — SM. The magnitude of the flow coefficients generally decreases with increas- 
ing n. There is a local maximum at r/ ^ 0. While the p+p results are approximately 
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symmetric, as expected, the p+Pb results exhibit a strong asymmetry. In the direc- 
tion of the proton beam, the p+Pb coefficients increase to approach the magnitudes 
of the p + p coefficients. 




Fig. 7. The flow coeflicicnts i>2{2} (left) and t>3{2} (right) as a function of ccntrality in p+Pb 
collisions at (r;) = 4, 0, and —4. The results are shown for both AMPT - def and AMPT - SM. 



Figure [7] shows the centrahty dependence of W2{2} and W3{2} in p+Pb coUisions. 
Three different 77 ranges are shown: —5 < rj < —3, \'q\ < 1, and 3 < i] < 5 {{rj) = —4, 
and 4 respectively. The coefficients increase from central to peripheral events for 
both AMPT - SM and AMPT - def. The coefficients are largest at (77) = while those 
at (77) = 4, in the direction of the Pb beam, are smallest, consistent with the results 
in Fig. [6l 

3.3. Medium Modification Factor RpPh 
3.3.1. RpPhiPx) at r] ^ 

We now turn to the medium modification or suppression factors, defined as the 
ratio of the p + A to p + p cross section, 

Rt-hjPryb)^ 1 dapPh/drjdpT ^^^^ 

Results are shown first for RpPhipr) for charged hadrons (pions) in Fig.[8l Standard 
shadowing calculation are shown on the left-hand side, results with event generators 
in the center, while CGC calculations and results with energy loss in cold matter 
are shown on the right-hand side. 



Cronin effect expected at the LHC: dipole formulation (B. Z. Kope- 
liovich and J. Nemchik) At LHC energies, hadron production is dominated by 
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Fig. 8. Charged particle RpPhPT at y/s~^^ = 5.02 TeV at ~ 0. (Top left) Result s with m ore 
'standard' shadowing (labeled EKS98, EPS08 and HKN) described in Sec. [23 Refs.EHIIMl (la- 
beled Kopeliovich) and the HIJING2.1 shadowing parameterization in Eq. II21I I are compared. The 
difference in the HIJING2.1 curves depends on whether the hard scatterings are coherent or not, 
see Sec. 12.31 (Top Right) HIJINGBB2.0 with and without shadowing (Sec. I2.4|l compared to AMPT 
default and with string melting (Sec. [231l . (Bottom) The band from rcBK saturation model cal- 
culations by Albacete et al. and Rezaeian with N = 5 and varying aj," described in Sec. 12.11 are 
compared to IP-Sat calculations by Tribedy and Venugopalan in Sec. l2.2l and calculations by Vitev 
et al. discussed in Sec. 12.61 (More detailed results for the uncertainties in Rezaeian's calculation 
at other rapidities can be found in Fig.[TOl) The ALICE results from Ref.l^l are also shown. The 
systematic uncertainties are shown in blue, the statistical uncertainties are in black. 



gluon fragmentation. The dipole formalism, on the light cone, is employed in the 
target rest frame, leading to the factorized expression for hadron production: a 
convolution of the projectile gluon PDF, the gluon fragmentation function, and the 
gluon splitting, g gg, cross section.'iSll The gluon splitting cross section can be 
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written in the dipole representation as 
d(j{gA ^ gigX) 



1 / ^Nrp (i\ 1 „N /- 



d'^PT dyi 



where the subscript 1 indicates the hadronizing gluon involved in inclusive high pT 
hadron production by gluon fragmentation. The momentum fraction of radiation 
gluon gi is defined as a = the transverse coordinates of the emitted gluons 

are ri and r2, x is the momentum fraction of the gluon in the target nucleus 
and TAib) is the nuclear thickness function. The distribution function, ^gg{r,a), 
derived in,QSSl describes the \gg) Fock component of the projectile gluon light-cone 
wavefunction, including the nonperturbative gluon interaction. It is characterized 
by the localization of gluons at a short relative transverse separation, ~ 0.3 fm. 
The ig dipole cross section, cr|^, is expressed as a combination of the qq dipole 

cross sections, extracted from phenomenology and DIS data, cr^gir) ~ fl^rggl^) + 

<yqqiar) + (7qq[{l - a)r]^ 

Multiple interactions of higher Fock components (containing more than one 
gluon) of the incident proton in the nucleus leads to a suppression of the dipole 
cross section, known as gluon shadowing. The magnitude of the suppression factor 
Rg was evaluated in Ref.'^^ and, as a consequence of the small gluon separation, tq, 
the modification was found to be rather small. Thus, at first order, gluon shadowing 
can be implemented by the simple replacement — >■ Rg{x, Q^, b)a-ig in Eq. ([M]). 
The results Ref.liSSJ a,t ^ — 5 TeV are given by the solid cyan curve on the left-hand 
side of Fig. [51 

The Cronin enhancement is rather weak with a maximum of RpPb ~ 1.05 at 
Pt ~ 3 GeV/c. The height of the maximum is extremely sensitive to the strength 
of gluon shadowing. The enhancement could completely disappear, or even change 
sign, if the strength of gluon shadowing was underestimated in Ref.^i^SI 

Results with standard shadowing parameterizations (G. G. Barnafoldi, 
J. Barette, M. Gyulassy, P. Levai, G. Papp and V. Topor Pop) The left- 
hand side of Fig. [8] also shows RpPhiPr) for several standard shadowing functions: 
EKS9#2l (blue), EPS08^3 (magenta), and HK^P (green). All three were calcu- 
lated in the rapidity interval \r]\ < 0.3, in the center of mass frame. Since none of 
these parameterizations include impact parameter dependence, the results are in- 
dependent of collision centrality. For details, see Ref.'^ For other results employing 
kTpQCD.v2.0, see Refs.lHUIMEnZEniEII] 

HIJING2.1 (R. Xu, W.-T. Deng and X.-N. Wang) Two HIJING2.1 calcula- 
tions including shadowing with parameter Sg — 0.28 are shown with the standard 

'The X dependence of has been suppressed. 
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shadowing parameterizations on the left-hand side of Fig. [8l They are lower than 
the other calculations at high px ■ Including the decoherent hard scatterings moves 
RpPh closer to unity over all px- 

HIJINGBB2.0 (G. G. Barnafoldi, J. Barette, M. Gyulassy, P. Levai, M. 
Petrovici and V. Topor Pop) The HIJINGBB2.0 results without (NS) and with 
(WS) shadowing in the central rapidity region, \ri\ < 0.8, are shown in the central 
panel of Fig. [8] Including shadowing and strong color field effects reduces RpPh 
from unity to ^ 0.7 for 5 < < 10 GeV/c. This is similar to CGC predictions in 
the KKT04 model. Further, similar, results for Rpphipr) are found in Ref.liiSlEil] 
using LO pQCD coUinear factorization with the HIJING2.0 shadowing parameteri- 
zation,!^ the GRV proton PDFs,'^and hadron fragmentation functions from Ref.l^ 

AMPT (Z. Lin) The AMPT results, evaluated in the center of mass frame pseudora- 
pidity interval |r/| < 1 are also shown in the central panel of Fig. [S] In minimum bias 
collisions, A^coii — 7.51 in AMPT. There is an w 50% suppression of charged hadron 
production for the entire pr range shown, considerably lower than that obtained 
in the other calculations shown. The statistical uncertainty, not shown in Fig. |8l 
becomes large for pT > 10 GeV/c. There is little difference in RpPh between the 
default and string melting options of AMPT. 

The right-hand side of Fig. [8] compares rcBK results from Albacete et al. and 
Rezaeian to IP-Sat calculations by Tribedy and Venugopalan and pQCD calcula- 
tions by Vitev et al. including energy loss in cold matter. 

rcBK predictions (J. Albacete, A. Dumitru, H. Fujii, Y. Nara, A. Reza- 
eian and R. Vogt) Rezaeian's calculation assumes iVcon = 6.9 in minimum bias 
p+Ph collisions .ISH To compare these calculations to data, it is necessary to rescale 
the results to match the normalization A^'cou to the experimental value. The NLO 
MSTW proton PDF^ISl and the NLO KKP fragmentation functions^^ are em- 
ployed. The rcBK equation, Eq. ([7]), accounts for the rapidity /energy evolution of 
the dipole but does not include impact parameter dependence explicitly. In the 
minimum-bias analysis here the initial saturation scale, Qos{xq = 0.01), is con- 
sidered as an impact-parameter averaged value and extracted from minimum-bias 
data. Thus, fluctuations in particle production and any other possible nonpertur- 
bative effects are incorporated into this average value of Qos- It was found that 
Qlp{xQ = 0.01) « 0.168 GeV^/c^ with 7 sa 1.119 provides a good description of 
small-x proton data from the LHC, HERA and RHIC,'!^^^ see also Ref.tS^ These 
values are employed in the rcBK description of the projectile proton. 

Here the difference between protons and nuclei originates from the different 
initial saturation scales, Qos, in the rcBK equation, see Eq. (|8]). The initial nuclear 
saturation scale Qqa at a given centrality is generally less constrained than that in 
the proton because the small-x data from d+Au collisions at RHIC and DIS data on 
heavy target are limited with rather large uncertainties, leading to correspondingly 
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larger uncertainties on Qqa^ 

(3 - 4) Ql^ < QlA^o - 0.01) < (6 - 7) Ql^ . (39) 

The role of fluctuations on particle production and other nonperturbative effect, 
including the impact-parameter dependece of the initial nuclear saturation scale 
are effectively incorporated into the average value of QoA assuming Qq^ = ^Qop 
is constrained by experimental data. Figure |S] shows results at 77 = with N = 5 
and varying a". For results over the full range of N in Eq. ([5^ . see Fig. (TU] 

On the other hand, Albacete et al let the initial nuclear saturation scale be 
proportional to the nuclear density at each point in the transverse plane. They 
compare two different methods,'^ 

QlAib,xo) = Np,,t{b)Qlp{xo) and (40) 

= Nli^,ib)Ql^{xo) . (41) 

The second option is an ad hoc way of correcting for the violation in the additivity of 
the nucleon number in the dipole scattering amplitude at small dipole sizes resulting 
from the fact that 7 7^ 1. The position of the nucleons in the transverse plane 
is simulated by Monte Carlo methods makes it possible to account for geometry 
fluctuations which can have large numerical impact.'^ 

IP-Sat (P. Tribedy and R. Venugopalan) The nuclear modification factors 
for inclusive charged particles in minimum-bias collisions employing the IP-Sat 
approach is done using fc^-factorization approach at y = in the center-of-mass 
frame. The LO KKP fragmentation function^^ are used to compute the inclusive 
charged hadron spectrum from the gluon distribution. The details of the parameters 
used can be found in Refs.'^'^ The band shows the uncertainty in the calculation. 
Note that Rpph approachs unity with increasing px- 

Cold matter energy loss (Z.-B. Kang, I. Vitev and H. Xing) The red 

curves on the right-hand side of Fig. [8] are results including cold matter energy loss 
at y = 0, see Sec. 12.61 The upper curve corresponds to parameters determined for 
RHIC while the lower curve includes a potential enhancement in these parameters 
for LHC energies. There is a very small "Cronin peak" in the low pT region. The 
peak is very close to unity and not as pronounced as that seen in lower energy 
fixed-target experiments because dynamical shadowing strongly suppresses particle 
production in this region. Initial-state energy loss is also larger because diagrams 
with incident gluons make a bigger contribution to the cross section. At high px, a 
small, ~ 15 %, suppression remains, also due to cold nuclear matter energy loss. 

The preliminary data are most compatible with the saturation model calcula- 
tions, albeit the uncertainties of these calculations are still rather large. The calcu- 
lation by Kopeliovich is also in relative agreement with the data. However, we note 
that none of the calculations available for both dNch/drj and Rpph at 77 ~ agree 
with both observables simultaneously. 
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3.3.2. HIJING2.1 parton vs. hadron RpA (R- Xu, W.-T. Deng and X.-N. 
Wang) 

On the left-hand side of Fig. [SI the nuclear modification factor of the final-state 
parton pT spectra, Rpa{pt), is shown for p+Ph collisions at the LHC. Here (Ncow) 
is the average number of binary nucleon-nucleon interactions in p + A collisions. It 
is clear that both nuclear shadowing and the hard-soft coupling can suppress the pr 
spectra of the produced partons. These features will transferred to the final hadron 
spectra after hadronization, as shown on the right-hand side of the Figure. 
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Fig. 9. (Left) The nuclear modification factor of the parton pT spectra in p-l-Pb coUisions. (Right) 
The charged hadron nuclear modification factor with different HIJING2.1 options. The arrow in- 
dicates the most probable trend of the nuclear modification factor to transition from the low to 
the high regions. 



The predicted charged hadron nuclear modification factors are shown on the 
right-hand side of Fig. |9l The enhancement at intermediate pt is due to the kx 
broadening arising from multiple scattering. Parton shadowing, the modified soft- 
hard coupling and enhanced density of gluon jets due to valence quark conservation 
all contribution to suppression of the charged hadron spectra. In addition, string 
fragmentation can further suppress the high pr hadron spectra compared to inde- 
pendent fragmentation. 

Since the effects of parton shadowing will disappear at large due to the 
QCD evolution not yet included here while hard and soft scatterings will decohere, 
the nuclear modification factor will likely follow the default result, including DHC, 
at low Pt- At large pT, the result should approach that with DHC effects but no 
shadowing. There may be further possible modifications due to the hadronization 
of multiple jets. In Fig. [3J this probable trend is indicated by the arrow joining 
these two results at intermediate pr- 

3.3.3. RpPhipT) at \rj\ ^ 

Here we show several results for the charged particle RpPb away from central rapid- 
ity. Two employ the rcBK approach but with different initial conditions. The last 
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is a pQCD calculation including cold matter energy loss. 



Rezaeian rcBK (A. Rezaeian) In Fig.[TUl R^{pt) is shown for minimum bias 
collisions at 77 = 0, 2, 4, and 6 (with the convention that the proton beam moves 
toward forward rapidity) obtained from the hybrid factorization Eq. (|H) supple- 
mented with rcBK evolution. The band labeled CGC-rcBK includes uncertainties 
associated with the variation of initial nuclear saturation scale, Qq^, constrained in 
Eq. (P^ . as well as uncertainties due to the choice of factorization scale, = 2pT, 
PT, and pr/2, in Eq. (|4]). The values of N scale the nuclear saturation scale relative 
to that in the proton, Qq^ = NQ^^, with 3 < < 7. The variation of the results 
with N represents an effective variation with centrality and incorporating possible 
effects of fluctuations. Going forward in proton rapidity reduces the range of the 
results and lessens the dependence on N. 

As shown in Ref.,'^the RHIC data unfortunately cannot fix the value of a!," 
in Eq. (jH): with a reasonable ii'-factor, < a!," < 0.3 is consistent with hadron 
production at RHIC both mp + p and d+A collisions. Figure [TO] shows the effect of 
varying a™ in Eq. Q. For each value of r/, results are shown both for a'^ = and 
variation within the range 0.09 > a!," > 0.3. For r/ = and 2, both a'/" = 0.1 and 
0.2 are shown. There is only a small difference at 77 = while a larger difference can 
be seen at forward proton rapidity, especially for N > 5. With rj = A and 6, only 
a™ = 0.1 is shown. Note that taking a™ > changes the pt slope considerably, as 
may be expected. 

In order to quantify the uncertainties due to different initial nuclear saturation 
scales, in Fig. [TO] is shown for different initial saturation scales Qg^ = NQ^p 
ioi 3 < N < 7 with factorization scale Q = pt- Unfortunately, is very sensi- 
tive to Qqj^- However, by comparing to measurements of at one rapidity, the 
predictions in Fig. [TO] can be used to determine Qqa that rapidity. Then, if the 
measured Rp\ at one rapidity lies between two lines labeled Ni and N2, the results 
at other rapidities should correspond to the same values of N and thus the same 
QIa- Only these results should be considered to be true CGC predictions. In this 
way, knowing R'^ at one rapidity significantly reduces theoretical uncertainties as- 
sociated with at other rapidities. Therefore, despite rather large theoretical 
uncertainties, it is still possible to test the main CGC/saturation dynamics at the 
LHC. 

It is clear that CGC approaches predict more suppression at forward proton ra- 
pidities compared to collinear factorization.liilEill Moreover, the small-x evolution 
washes away the Cronin-type peak at low pT at all rapidities. Therefore, observa- 
tion of a Cronin peak in p+Ph collisions at the LHC, regardless of whether it is 
accompanied by overall enhancement or suppression, can potentially be a signal of 
non-CGC physics since it is difficult to accommodate this feature in the current 
accuracy of the CGC approach.lSSl 
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Fig. 10. The nuclear modification factor for charged hadron production in minimum bias 

p+Pb collisions at ?7 = 0, 2, 4, and 6 (with the convention that the proton beam moves toward 
forward rapidity) obtained from the hybrid factorization in Eq. I I39I1 assuming different values of 
the saturation scale in the nucleus, Qq^- The lines labeled by a given value of Af, for 3 < A'^ < 7, 
are results with fixed factorization scale p^p = pT and fixed saturation scale Qq^ = NQ'^^ and 
Qqp = 0. 168 GeV^ / . The bands shown the variation in the results with the choice of factorization 
scale. Two panels are shown for each rapidity. The upper panel shows results obtained by taking 
Cfg" = (assuming only elastic contribution) while the bottom panel shows the variation of in 
the range 0.09 > a™ > 0.3. In the bottom panels for i) = and 2, results with both a^^ = 0.1 and 
0.2 are shown, while for rj = i and 6, only a^" = 0.1 is shown. The plots are taken from Ref.l^ 
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Albacete et al rcBK (J. Albacete, A. Dumitru, H. Fujii and Y. Nara) 
Figure [TT] shows bands for Rpa{pt) at r] — and 2 with rcBK-MC. The upper 
edges of the bands are calculated with 7 = 1.119 and KKP-LO fragmentation 
functions, the lower limit of the bands corresponds to the McLerran-Venugopalan 
initial condition (7 = 1) with DSS-NLO fragmentation functions, taken from Ref.l^ 
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Fig. 11. The nuclear modification factor for three different centrality classes assuming fcy- 
factorization. The 77 = 2 result is obtained with the convention that the proton beam moves 
toward forward rapidity. 



There is suppression of RpPh{PT IGeV/c) = 0.6 ± 0.1 at midrapidity. The 
ratio increases monotonically to pt ~ 3 GeV/c where it flattens. There is no clear 
Cronin peak. For all unintegrated gluon densities, RpPh decreases with increasing 
rapidity. In addition to the minimum bias results, two separate centrality classes 
are also shown, a central bin, A^part > 10, and the most peripheral bin, iVpart < 5. 
There is a stronger effect for the most central collisions, as expected, for low pT 
at both rj — and 2. (Note that there is a smaller difference between minimum 
bias and the most central collisions at 77 = 2. The overall effect is very weak for 
the peripheral bin with suppression persisting only to pr 2 — 3 GeV/c. For more 
details see R.ef 

Cola inatter energy loss (Z.-B. Kang, I. Vitev and H. Xing) Figure [12] 
presents model predictions for RpPhiPx) in minimum bias p+Ph collisions at y = 
(top), y — 2 (middle), and y = 4 (bottom). The upper edge of the bands corresponds 
to the RHIC scattering parameters while The lower edge allows for a potential high- 
energy enhancement of the parameters. 

The y = results were already presented in Fig. [H At large proton rapidity 
the CNM effects are all amplified due to the larger values of the proton momentum 
fraction X2 (relevent for cold nuclear matter energy loss with the proton moving 
in the direction of forward rapidity), the smaller values of the nuclear momentum 
fraction xi (relevent for dynamical shadowing with the nucleus moving in the di- 
rection of backward rapidity) and more steeply-falling px spectra (relevent for the 
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Fig. 12. Predictions for the nuclear modification factor Rppb as a function of pT for charged 
hadron production in minimum bias p+Pb collisions. Results are shown for three rapidities: y = 
(top), y = 2 (center), and J/ = 4 (bottom) with the convention that the proton beam moves toward 
forward rapidity, see 

Cronin effect). As a result, at low px and large ion rapidity, dynamical shadow- 
ing can dominate and lead to stronger suppression of inclusive particle production 
(note the disappearance of the small Cronin enhancement at y > 0). At high pr the 
suppression is a combined effect of cold nuclear matter energy loss and the Cronin 
effect. As will be seen later, the behavior of RpPh for tt'~' and direct photon pro- 
duction is qualitatively similar to the RpPb dependence for charged hadrons shown 
here. 

Away from midrapidity, these predictions suggest a stronger effect than that 
found with EPS09 parton shadowing alone. On the other hand, the effect shown 
here is weaker than the results using the rcBK-MC approach also shown in this 
section. 



AMPT at J7 ^ (Z. Lin) Figure [T3l shows results with the AMPT — def calculation in 
minimum bias p+Ph collisions at five different r/ values. At r/ = 4, the ratio is larger 
than at lower rapidities. The ratios are very similar for rj = and 2. The backward 
rapidity results have lower values of i?ppb- While the ratios are indepedent of px 
within the uncertainties for 77 = — 2, and 2, the ratios increase with px for r/ = 4 
and —4 with the largest increase for rj — —4. 
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Fig. 13. Predictions for the nuclear modification factor Rppb as a function of pT for charged 
hadron production in minimum bias p+Pb coUisions for rj = —4, —2, 0, 2 and 4 calculated using 
AMPT - def . 



3.3.4. Forward-backward asymmetry (G. G. Barnafoldi, J.Barette, M. 

Gyulassy, P. Levai, G. Papp, M. Petrovici, V. Topor Pop. Z. Lin and 
R. Vogt) 

One method of determining the difference between the proton and lead rapidity 
regions is studying the forward-backward asymmetry in charged hadron production, 
given by 



'd'^PTdri\rj<o 
I > 0^ 

(42) 



Ehd^a'^pJd^PTdr^l 
^pPb(PT,r/>0) 



^pPb(PT,r/<0) ■ 

where the 'forward' direction, rj > 0, is that of the lead beam, toward positive 
rapidity while the 'backward' direction, 77 < is that of the proton beam, toward 
negative rapidity. 

The asymmetries, calculated in the center of mass frame in the range 0.3 < jryj < 
0.8,1^^ both for standard shadowing and event generators, are shown in Fig. [141 Two 
calculations of the forward-backward asymmetry in p-|-Pb collisions are shown. 
The first result is in the kTpQCD_v2.0 approach with several standard shadowing 
parameterizations. The second is obtained with AMPT. 

Results with kTpQCD_v2.0 including the HIJING2.0,''^ EKSDS,"^ EPS08,'22l and 
HKN^'* shadowing parameterizations are presented. Since the EKS98, EPS08 and 
HKN parameterizations are independent of impact parameter, there is no difference 
between minimum bias results and results with a centrality cut. Thus, these are 
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all labeled MB for minimum bias. With the HKN parameterization, Yasym ~ 1 
independent of pt- The EKS98 and EPS08 parameterizations, producing identical 
ratios, predict a small enhancement, Yasym > 1- 
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Fig. 14. Predictions for the forward- backward asymmetry, V^sym (pt), from Refs.E^EHlcentrality 
independent results are shown for the HKN, EKS98 and EPS08 parameterizations (labeled MB). 
Minimum bias results are also shown for HIJINGBB2.0 and HIJING2.0 with multiple scattering. In 
addition, HIJING2.0 results in MB collisions and for the 20% most central collisions (using the 
parameter Sa in Eq. I I22II ) are also shown. The blue points are the AMPT — def results. 



Results with the HIJING2.0 shadowing parameterization with and without mul- 
tiple scattering differ for pT < 5 GeV/c but converge at higher pr- When the 
collision centrality is taken into account, the forward-backward asymmetry is en- 
hanced by the centrality cut alone. Chaning the strength of low x shadowing, as in 
Eq. ([22|) . does not change the shape of Yasym (pt) for pr > 5 GeV/c. 

Finally, results for the minimum bias asymmetry is also shown for 
HIJINGBB2.0.'^ These calculations include both shadowing and strong color field 
effects (indicated SCF"). Here the ratio rises from near unity at ^ to a 
peak at pr ~ 5.5 GeV/c. It then decreases to Yasym ^ 0.75 for pr > 10 GeV/c. 

The forward-backward asymmetry of i?ppb, calculated in the center of mass 
frame with AMPT — def in the range 0.3 < \ri\ < 0.8, is shown by the blue points in 
Fig. [M] for charged particles in minimum-bias events with simulations employing 
AMPT — def. The initial enhancement first increases with px and then decreases. 
Although the statistical errors become very large above px ^ W GeV/c, it is clear 
that the result is similar to the general behavior of the HIJINGBB2.0 result. 

The low Pt enhancement may arise because the asymmetry in dN/ drj in p+Pb 
relative to p + p in the lead direction (positive rapidity) introduces a natural en- 
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hancement in i?ppb in the forward direction. The enhanced particle production in 
p+Pb is most hkely at low pr, see e.g. Fig. [T3l 

The AMPT — def results include the impact-parameter dependent nuclear shad- 
owing parameterization implemented in HIJINGl.O. Therefore the asymmetry at 
other centralities would differ from that in minimum-bias events. 

3.3.5. RpPbirj) (G. G. Barnafoldi, J. Barette, M. Gyulassy, P. Levai, M. 
Petrovici and V. Topor Pop) 

Finally, we end this section with the predicted ratio iZpPb as a function of pseu- 
dorapidity in HIJINGBB2.0. As discussed previously, increasing the minijet cutoff 
parameter po from 2 GeV/c va. p -\- p collisions to 3.1 GeV/c in p+Vh causes a 
slow growth in the per nucleon multiplicity which could be interpreted as evidence 
for gluon saturation. It is difficult to directly relate po to the saturation scale Qs 
because in HIJING hadronization proceeds through longitudinal string fragmenta- 
tion. The low pt part of the minijet spectrum is particularly sensitive to the ^s„„ 
and A dependence of minijet suppression while the > 5 GeV/c minijet tails are 
unaffected. The energy and A dependence of the string tension arises from strong 
color field (color rope) effects not included in CGC phenomenology that assumes 
hadronization by fcr-factorized gluon fusion. 
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Fig. 15. The suppression factor -Rppb as a function of pseudorapidity, ri, for HIJINGBB2.0 with 
strong color field effects, without (solid) and with (dashed) shadowing. 

HIJING minijet hadronization does not proceed by independent fragmentation, 
as in PYTHIA,'SSl but through string fragmentation with gluon minijets represented 
as kinks in the strings. The interplay between longitudinal string fragmentation 
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dynamics and minijets is a nonperturbative feature of HIJING-type models. The 
effect of string fragmentation on the multipHcity is manifested in the behavior of 
RpPhiv) the fragmentation regions |r/| > 5. In the nuclear beam fragmentation 
region, the ratio is approximately linearly increasing from ry = 5 to the kinematic 
limit. There is an approximate plateau in the ratio over < ry < 5, followed by a 
decrease toward the value l/A'^coii in the proton fragmentation region, a Glauber 
geometric effect first explained in Refs-l^i^lii^ and a feature of string fragmentation 
in HIJING. 

4. Identified Particles 
4.1. AMPT (Z. Lin) 

Some representative AMPT results for identified particles are shown here. The ra- 
pidity distributions for if + and production in the default and string melting 
versions are shown in Fig. 1161 The default and string melting versions of the rapidity 
distributions of protons and antiprotons are shown in Fig. 1171 

There are differences between the default and string melting versions, even in 
p + p collisions. However, these differences are typically small. Note that the bump 
at midrapidity in the proton and antiproton rapidity distributions in the string- 
melting version is likely a result of the simple quark coalescence model implemented 
in AMPT (see e.g. Fig. 38 of ReiW^. 




y y 
Fig. 16. Rapidity distribution, dN/dy, of (left) and (right) mesons. 

The left-hand side of Fig. [18] shows RpPh for protons, antiprotons, neutral pions 
and charged particles. The identified particles are shown for y — while the charged 
particle ratio is for 77 = 0. The ratios are all very similar except for protons which 
grows almost linearly with pT- 

The right-hand side of Fig. [18] presents results for the suppression factor for 
charged particles within \ri\ < 1 in three centrality bins: minimum bias; the 10% 
most central collisions; and a peripheral bin, the 60 — 80% most central collisions. 
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Fig. 17. Rapidity distribution, dN/dy, of and p and p in the default (left) and SM (right) scenarios. 



Interestingly, the ratio is smallest and decreasing with pt for the more peripheral 
bin while the largest value of RpPb is the most central bin. 




Fig. 18. (Left) The minimum bias results at y = for tt", p and p arc shown with the charged 
hadron ratio for rj = 0. All results are calculated with AMPT — def . (Right) The charged hadron 
results at midrapidity, |r;| < 1, are shown for minimum bias collisions as well as in two centrality 
bins: the 10% most central and the 60-80% most central. All results are calculated employing 
AMPT -def. 



4.2. Centrality dependent Rppi,{PT,y = 0) with EPSOQs nPDFs (I. 
Helenius and K. J. Eskola) 

Predictions of the centrality dependence of tt^ production at midrapidty are pre- 
sented here. Two spatially-dependent nuclear PDF (nPDF) sets, EPS09s and 
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EKS98s were determined in RefP^ The key component is a power series ansatz 
of the nuclear thickness function (Tj^^s)) for the spatial dependence of the nPDF 
modifications, 



' 1 



(x, Q',s} + Q') i^AisjY . (43) 



The A dependence of the earlier global nPDF fits EPSO#Sl and EKS9S'^ was ex- 
ploited to obtain the values of the A- independent coefiicients Cj{x, Q^). It was found 
that n = 4 is sufficient for reproducing the A systematics of the globally-fitted 
nPDFs. 

Predictions are presented for the nuclear modification factor -Rpp^ of inclusive 
neutral pion production at the LHC using the new EPS09s NLO nPDF set. As in 
Ref . the centrality classes are defined in terms of impact parameter intervals. 
For a given centrality class, -Rpp^ is defined as 

\ dprdy / a^b / 

R^^^{pr,y;hM) = = -— ^, (44) 



cr-n dprdy •'"^ dprdy 



where the impact parameter limits hi and 62 are calculated from the optical Glauber 
model. The proton is assumed to be point-like, i.e. TpPb(b) = Tpb(6). The minimum 
bias ratio i^pp^ is obtained with bi = and 62 — >■ 00. The cross sections are 
calculated in NLO using the INCNLO packagdJ* with the CTEQ6M proton PDEspl 
along with three different fragmentation functions: KKpSa AKKfEIl and fDSS.Eni 
The uncertainty band reflecting the nPDF uncertainties is calculated using the error 
sets of EPSOQs with the fDSS fragmentation functions according to the prescription 
described in the original EPS09 paper.Esl 

In Fig. [19] the minimum bias Rppy^ is shown as a function of pt at midrapidity 
{y — 0) in the nucleon-nucleon center-of-mass frame. Both logarithmic and linear 
Pt scales are shown. In the region px < 10 GeV/c a suppression due to small- a; 
shadowing in the nPDFs is observed while at 10 < < 200 GeV/c there is a small 
enchancement due to the antishadowing. Even though the different fragmentation 
functions may yield different absolute cross sections, these differences cancel in the 
ratio i?ppb- 

In Fig. [20] -Rppb is presented in four different centrality classes, (0 — 20) %, 
(20 - 40) %, (40 - 60) % and (60 - 80) %. The impact parameter limits for these 
centrality classes can be found inRe^Ii^niFor comparison, the minimum bias results 
are also shown. Slightly larger nuclear effects are observed in the most central 
collisions than in the minimum bias collisions. Nuclear effects in the most peripheral 
collisions are about a factor of two smaller than those in the most central collisions. 



"http://lapth.in2p3.fr/PHOX_FAMILY/readmc_inc.html 
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Fig. 19. The nuclear modification factor for inclusive n production in p+Pb collisions at j/ = 
for minimum bias collisions calculated with the EPS09s NLO nPDFs, plotted on logarithmic (left) 
and linear (right) px scales. The blue error band is calculated employing the EPSOQs error sets. 
All scales are fixed to the pion pj-. Based on Ref.ESHl 




Fig. 20. The nuclear modification factor for inclusive n" production in p+Pb collisions at j/ = in 
four centrality classes calculated with the EPSOQs NLO nPDFs. The blue error band is calculated 
employing the EPSOQs error sets. All scales are fixed to the pion pj-. For comparison, the minimum 
bias result is also shown. Based on Ref.E^ 



4.3. Cold matter effects with energy loss (Z.-B.Kang, I. Vitev and 
H. Xing) 

The predictions for -Rpp^,, like those for charged hadrons shown in Fig. [Tulare based 
on perturbative QCD factorization. Cold nuclear matter effects are implemented 
separately. The advantage of this approach is that all CNM effects have a clear 
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physical origin, generally centered around multiple parton scattering. The imple- 
mentation of these calculable CNM effects is well documented, see Sec. 12.61 Isospin 
effects, the Cronin effect, cold nuclear matter energy loss, and dynamical shadowing 
are all included. 

The neutral pion results in Fig. [21] are rather similar to those shown in Fig. [121 
The upper edges of the bands at ?; = (top), y = 2 (middle), and y — 4: (bottom) 
correspond to the RHIC scattering parameters while the lower edges correspond to 
potential enhancement of these parameters. 
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Fig. 21. Predictions for the nuclear modification factor Rpp^^iPr) as a function in minimum bias 
collisions at y/s~^^ = 5 TeV. Results for three rapidities: y = (top), y = 2 (middle), and y = 4 
(bottom) are shown. The calculations have been made with the convention that the proton beam 
moves toward forward rapidity, see Ref.E3 



As seen by comparison with the EPS09 minimum bias results for y = in 
Fig. 1191 the combined effects included here are stronger than with shadowing alone. 
Antishadowing in EPS09 produces a ratio larger than unity for 10 < px < 200 
GeV/c while the ratio is smaller than unity over the entire range calculated here 
{pT < 50 GeV/c). 

5. Quarkonium (R. Vogt) 

The predictions for the J/ip suppression factor, considering only shadowing effects 
on the parton densities are described in this section. There are a number of pos- 
sible cold matter effects on J/iJ; production, including modifications of the parton 
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densities in nuclei (shadowing); breakup of the quarkonium state due to inelastic 
interactions with nucleons (absorption); and energy loss in cold matter. Since the 
quarkonium absorption cross section decreases with center-of-mass energy, we can 
expect that shadowing is the most important cold matter effect at midrapidity, see 
Refs.^^^ ^^'^ Here we show results for the rapidity and pt dependence of shadowing 
at ^/s^ — 200 GeV for d+Au collisions at RHIC and the rapidity dependence at 
^ySfjj^ = 5 TeV p+Ph collisions, neglecting absorption. 

The results are obtained in the color evaporation model (GEM) at next-to- 
leading order in the total cross section. In the GEM, the quarkonium production 
cross section is some fraction, Fc , of all QQ pairs below the HH threshold where 
H is the lowest mass heavy- flavor hadron, 

ag™(s) =Fc^ / ds / dxidx2 /f(xi,/.|) f^{x2,til) a.j{rs,fil,^il) ,(45) 

■ ■ J im^ J 

where ij = qq or gg and 'Jy(s) is the ij QQ subprocess cross section. The 
normalization factor Fc is fit to the forward (integrated over xp > 0) J/ip cross 
section data on only p, Be, Li, G, and Si targets. In this way, uncertainties due to 
ignoring any cold nuclear matter effects which are on the order of a few percent in 
light targets are avoided. The fits are restricted to the forward cross sections only. 

The same values of the central charm quark mass and scale parameters are em- 
ployed as those found for open charm, m = 1.27±0.09 GeV/c^, HF/m = 2.10lg'|5, 
and fiR/m = 1.60lg ijj-^^'' The normalization Fc is obtained for the central set, 
(to, fip/m, fiii/m) — (1.27 GeV, 2.1, 1.6). The calculations for the extent of the mass 
and scale uncertainties are multiplied by the same value of Fc to obtain the ex- 
tent of the J/ tp uncertainty band.li^il The results shown here are not the same as 
those calculated at leading ordeJ^^ because the LO and NLO gluon shadowing 
parameterizations differ significantly at low x^^ 

Figure [22] shows the uncertainty in the shadowing effect due to uncertainties in 
the EPS09 shadowing parameterizatiorpS (red) as well as those due to the mass 
and scale uncertainties obtained in the fit to the total charm cross section (blue) 
calculated with the EPS09 central set. All the calculations are NLO in the total 
cross section and assume that the intrinsic kx broadening is the same in p -I- p as 
in p-|-Pb. Note that the rapidity-dependent ratios are adjusted so that the lead 
nucleus moves toward negative rapidity and the small x region of the nucleus is 
at large negative rapidity, as is the case for pseudorapidity distributions and the 
rapidity-dependent ratios shown previously. 

The mass and scale uncertainties are calculated based on results using the one 
standard deviation uncertainties on the quark mass and scale parameters. If the 
central, upper and lower limits oi p /m are denoted as C, H, and L respectively, 
then the seven sets corresponding to the scale uncertainty are {(/zf/to, /^f/"^)} = 
{(C,C), iH,H), {L,L), {C,L), {L,C), {C,H), (i?,C)}. The uncertainty band can 
be obtained for the best fit sets by adding the uncertainties from the mass and scale 
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Fig. 22. The ratio -RdAu (left) and -RauAu (right) at ^/s = 5 TeV. The dashed red histogram shows 
the EPS09 uncertainties while the dot-dashed blue histogram shows the dependence on mass and 
scale. The p + p denominator is also calculated at 5 TeV and does not include any rapidity shift 
in p+Pb collisio ns. F or a discussion about normalizing the results to p + p collisions at different 
energies, see RefII2£ 



variations in quadrature. The envelope containing the resulting curves, 

^max — ^ccnt ^~ ^(*^/i,max *-^cciit)'^ ('^m.max '^ccnt)'^ ; (^^) 
f^min — f^ccnt \J (f^/A,min f^ccnt)^ ^" (^?n,min f^ccnt)'^ , (^'^) 

defines the uncertainty. The EPS09 band is obtained by calculating the deviations 
from the central value for the 15 parameter variations on either side of the central 
set and adding them in quadrature. With the new uncertainties on the charm cross 
section, the band obtained with the mass and scale variation is narrower than that 
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with the EPS09 variations. 
6. Photons 

Photons are ideal probes in heavy ion coUisions due to their lack of final-state 
interactions with either a quark-gluon plasma or a hot hadron gas. Baseline mea- 
surements in p + p and p+Ph collisions are very important to subtract photons 
from hard initial parton scatterings which are unrelated to QGP and, in particular, 
to determine cold nuclear matter effects on photon production. Measured nuclear 
modification factors for photons at high transverse momenta {pr > 5 GeV/c) mea- 
sured in Pb+Pb collisions at the LHC and Au-f Au collisions at RHIC, Raa ~ 1, 
have demonstrated that hard processes scale with the number of binary nucleon- 
nucleon collisions in nucleus-nucleus collisions. An especially important question 
for heavy-ion physics at the LHC is the magnitude of low-x effects on the parton 
densities, particularly for gluons. The presence of shadowing, i.e. suppression of 
the number of \ow-x partons, has a significant effect on benchmarking hard pro- 
cesses at intermediate px- Direct photons and other electromagnetic probes are in 
particular needed for a quantitatively precise determination of the nuclear parton 
distributions. 

This section includes an NLO calculation of direct photon production in p + p 
and p+Ph collisions; LO calculations of cold matter effects on photon production 
at several rapidities; saturation effects on direct and inclusive prompt photon pro- 
duction and photon-hadron correlations; and a discussion on gluon saturation and 
shadowing in dilepton and photon production. 

6.1. Direct photon cross sections (R. J. Fries and S. De) 

Predictions of the pT and y dependence of direct photon production are given 
here. The impact-parameter averaged EPS09 shadowing parameterization^^! ^j-g 
used with the CTEQ6.6 parton densitie^^ to calculate results in minimum bias 
p+Ph collisions. These calculations can thus be employed to check the validity of 
the EPS09 modifications. The per nucleon cross sections given here can be directly 
compared to the minimum bias p+Ph data by converting the cross sections to 
differential yields scaling the results by (Tin = 67 mb and the average number of 
NN collisions in a given centrality bin. 

These NLO in as calculations are performed with JETPHOXl. 3.1.2^2^ Hard 
prompt photons and fragmentation photons (obtained with the BFG-II fragmen- 
tation functionJi^ are both included. An isolation cut of i?T,hadion < 5 GeV on 
hadronic energy is imposed within a i? = 0.4 isolation cone around the photon, 
similar to the CMS analysis method ioi p + p collisions at 7 TeV. 

Figure [23] shows the rapidity dependence, da/dy, of direct photon production 
in p + p and p+Ph collisions, normalized per nucleon. A pT cut, px > 4: GeV/c is 
imposed. The results are presented in the lab frame. Recall that there is a difference. 
Ay = 0.465, between the center-of-mass and lab frames at the current energy. 
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Fig. 23. Tile direct plioton cross section, da/dy, as a function of rapidity in the laboratory frame 
for photons with pj- > 4 GeV/c and an isolation cut Erp hadron < 5 GeV. Both p + p and p+Pb 
results are shown. The p+Pb results are normalized to the per nucleon cross section. 

On the left-hand side of Fig.[24l the pr spectrum at y = 0, da/(PpTdy, is shown 
ioi p + p and p+Ph colhsions. The p + p result is scaled down by a factor of 100 
to separate the two curves. The corresponding nuclear modification factor is shown 
on the right-hand side. Note that a logarithmic scale is used on the a;-axis here to 
highlight the modification of the low pT part of the photon spectrum. 
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Fig. 24. (Left) The direct photon pT distribution at y = in the lab frame. The p+p distribution is 
scaled down by two orders of magnitude. (Right) The corresponding modification factor RpPbiPr)- 
Note the logarithmic scale. 
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6.2. Cold matter effects on photon production (Z.-B. Kang, I. 
Vitev and H. Xing) 

Prompt photon production in p + p collisions has two components, the direct and 
fragmentation contributions;®^ 



da _ dadir ^ da^ag 
dyd^PT dyd^pT dyd^pT 



The fragmentation contribution is given by: 



da. 



dyd'^P' 



frag ^K^Y.j ^d'kT,faMxi,k^Tj J ^ d' kr, hM^^, k^T.) 



dz 

X I ^D^f^{zc)Hat^c{s,t,u)S{s + t + u). (49) 



The expression is exactly the same as Eq. (1231) if the parton-to-hadron fragmentation 
function, Dfi/^{zc), is replaced by the parton-to-photon fragmentation function, 
D^/c{zc)- The direct contribution can be written as: 



dyd'^pi 



= [ !^(fkTja/N{xi,k'^ ) [ —d^kT2fb/N{x2,k'^ ) 

S J Xl J X2 

a,o 

xHab~^'y{s,i,u)S {s + i + u) , (50) 



where Hat^-y are the partonic hard-scattering functions for direct photon produc- 
tionPESni 

In p + p collisions, {k'^)pp — 1.8 GeV^/c^, along with a K factor of 0(1), which 
gives a good description of production at RHIC and LHC energies. The CTEQ6L1 
ppp^m are employed with the GRV parametrization for parton-to-photon frag- 
mentation functions.'^^ The factorization and renormalization scales are fixed to 
the transverse momentum of the produced photon, = /i/j = px- 

The results shown here are calculated employing the same cold matter effects 
described in Sec. 12.61 Figure [25] presents predictions for the nuclear modification 
factor in prompt photon production as a function of px for y ^ (top), y = 2 
(middle), and y = 4 (bottom). The upper edge of the bands corresponds to the 
RHIC scattering parameters The lower edge represents a high-energy enhancement 
of the parameters. The behavior of RpPb for direct photons is qualitatively the same 
as the TT^s and charged hadron results shown earlier. 
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Fig. 25. Predictions of the nuclear modification factor Hppb as a function of pT for prompt photon 
production in minimum bias p+Pb collisions for J/ = (top), y = 2 (middle), and 2/ = 4 (bottom), 
see Ref.113 



6.3. rcBK calculation of photon production (A. Rezaeian) 

The cross section for semi-inclusive prompt photon-quark production in p + A col- 
hsions at leading twist in the CGC formalism isjUHEHl 



p 



(Pbr (PpV dHr df]^ drjh V2{4n*) (pD'^V^ [p- h ~ l^-pV? 



y.5 



7 



xNFilh+PT^lxg), 



[21- p- It-Pt^ +p-{k- -p-)l'^ + l-{k- -l-){p^)^^ 



(51) 



where p"^ , I, and k are the 4-momenta of the produced prompt photon, the outgoing 
q (q) and the incident q (q), respectively. A K factor was introduced to absorb 
higher-order corrections. The light-cone fraction Xq is the ratio of the incoming 



quark and proton energies, Xq 



's/2. The pseudorapidities of the outgoing 



prompt photon, 77^, and quark, t]^, are defined as p~ = (pj./-\/2)e''^ and — 
{It / V2)e^'^ . The angle between the final-state quark and the prompt photon, A^, 
is defined as cos(A0) = {It ■ Pt) / {I'Tp't) ■ Only high pT light hadron production is 
considered here. Therefore rapidity and pseudorapidity are equivalent. The semi- 
inclusive photon-hadron production cross section in proton-nucleus collisions can 
be obtained by convoluting the partonic cross section, Eq. (|51|) . with the quark 
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distribution functions of the proton and the quark-hadron fragmentation function, 



dZf f , , , ^^qA^q(l)^{p-')X 



/ 1 



a ot d'^pT cI^lt drjj drjh 



where is the transverse momentum of the produced hadron. The sum over quark 
and antiquark flavors in Eq. (j52p is implied. The hght-cone momentum fractions 
Xq, Xq, and Xg in Eqs. (j5ip and (j52p are related to the transverse momenta and 
rapidities of the produced hadron and prompt photonj^^^il 

Xq =xq^^ (pI e"- + § e"^) , X, = ^ (p^ e^"^ + § e"' 

- Pll „min _ Pt ( e^h \ 

~ It ' ^/ ~ \ i~(p^/y/^)e''-' J 

The single inclusive prompt photon cross section in the CGC framework can be 
obtained from Eq. (1511) by integrating over the momenta of the final state quark or 
antiquark. The corresponding cross section can be then divided into two contribu- 
tions: fragmentation (first term) and direct (second term) photons ,1^^ 

,9^7, = / dXqfq[Xq,Q^)-NF{Xg,plp/z)D^lq{z,Q^) 

d-^brd'^PT drj-f I^tt; Lj^,„i„ z 

+ /^"'" / dXqfq{Xq,Q'')z''[l + {l-zf] [ d^T l^T {x g , h) 



2^2 (p; 



li,<C 



where D^/g{z ,Q^) is the leading order quark-photon fragmentation func- 
tion.l^'i^SESI] Similar to the hybrid formalism for inclusive hadron production, 
Eq. Q, Q is a hard scale. Although the cross sections given in Eq. ^ and Eq. (|53l) 
describe different final-state particle production, they are strikingly similar. The 
light-cone momentum fractions Xg, Xg, and z above are related to the transverse 
momentum and rapidity of the prompt photon,'^^ 



hg ~ XqC , U,g ~ ^ ^ 



Ml! _L (iT-pif 

z 1-z 



mill _ Pr^pTi-i r, — Pt ^ri^ 

The expression in Eq. ([55]) was obtained using a hard cutoff to subtract the collinear 
smg 

ularity.liSil The use of a cutoff may result in a mismatch between the finite cor- 
rections to Eq. (j53p and those included in parameterizations of the photon fragmen- 
tation function. However, this mismatch is higher-order in the coupling constant 
and its proper treatment requires a full NLO calculation. 

In Eqs. ((52|) and ([SS)) . the factorization scale /xf is assumed to be the same 
in the fragmentation functions and the parton densities. In order to investigate 
the uncertainties associated with choice of /ij?, several values oi fip are considered: 
fiF — 2pJ; pj; and p^/2. The amplitude Np in Eqs. (|5T|) and (f53| is defined in 
Eq. ©. 
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Fig. 26. The nuclear modification factor for direct photon production in minimum bias p+Pb 
collisions at rapidities rjj = 0, 2, 4, and 6 (with the convention that the proton beam moves toward 
forward rapidity) obtained from Eq. I I53I1 with solutions of the rcBK equation with different initial 
nuclear saturation scales. The band labeled CGC-rcBK includes uncertainties due to the variation 
of the nuclear saturation scale and the factorization scale fip . Similar to Fig. 1101 the lines labeled 
N are results with a fixed factorization scale, fip = pip, and a fixed saturation scale Qg^ = N'Qop 
with QL = 0.168 GeV^/c^. The range 3 < iV < 7 is constrained in Eq. l|39ll. See Ref.E^ 



In Figs, inland [271 predictions of the direct photon and single inclusive prompt 
photon production nuclear modification factors Rj^ijiT) in minimum bias p+Ph 
collisions at several rapidities. The solutions of the rcBK equation were obtained 
using Eq. (j53p with different initial nuclear saturation scales. The band labeled 
CGC-rcBK includes uncertainties due to the variation of Qg^ within the range 
given in Eq. with factorization scales iip — 2p7p, p7p, and p'^/2. In Fig. [271 the 

= results include inclusive prompt photon production calculated in Ref.li^sj 
employing the lancu-Itakura-Munier (IIM) saturation model which also provides a 
good description of the HERA data.l^^In the IIM method, saturation is approached 
from the BFKL region. Therefore the IIM small- a; evolution is different from that 
obtained with the rcBK equation. 

As discussed in Sec. 13.3.31 the results in Figs. [26l and [271 on -R^^, together with 
those on i?^^ in Fig. [TOl can be used to fix the nuclear saturation scale. Once Qqa 
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Fig. 27. The nuclear modification factor for single inclusive prompt photon production in 

minimum bias p+Pb collisions at ?7-y = and 2 (with the convention that the proton beam 
moves toward forward rapidity). The curves are described in Fig. 1261 T he d ashed red line labeled 
CGC-IIM was calculated according to the IIM dipole saturation model.LiS^I See Ref.[ini 

has been established at one rapidity, the predictions at other rapidities are fixed 
and the CGC/saturation dynamics can be tested. 



6.3.1. Photon-hadron correlations 

Photon-hadron aziniuthal correlations in p + A and p + p collisions could be an 
excellent probe of small- a; dynamics IIMIIISII -pji jg correlation can be defined aSj^i^^^^ 



P(A0) 



(Phx pl^dp!^ pTpdpTp drjj drjh i 



d^bx p!^dp!^ pTpdpTj^ dr]^ drjh d(j) 



.(54) 



The correlation function P(A0) is the probability of semi-inclusive photon-hadron 
pair production in a certain kinematic region at angle relative to production in 
the same kinematics at fixed reference angle, A0c = 7r/2.U^I^^ 

Fig. [25] shows the predicted P{A(p) in minimum bias p + p and p-|-Pb collisions 
at 5 TeV for pt^ < p7p (left) and > p7p (right). Given the pT of the produced 
photon and hadron, the corresponding correlation can have either a double or single 
peak structure. In Ref.^^^ it was shown that this feature is related to saturation 
physics and is governed by the ratio p^/p^. The change from a double to a single 
peak correlation, depending on the relative px, is unique to semi- inclusive photon- 
hadron production: since the trigger particle in dihadron correlations is always a 
hadron, it consistently exhibits a single-peak structure. 

Photon-hadron correlations can also be quantified by the coincidence probabil- 
ity. In contrast to production of a more symmetric final state such as dihadron 
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production, in photon-hadron production the trigger particle can be either the 
prompt photon or the hadron.E^H When the photon is used as the trigger, the coin- 
cidence probabihty is defined as CPh{^(f>) = Nf^'^" {A(f)) / where iV^''"^(A0) is the 
photon-hadron yield. The momentum of the photon trigger, the leading (L) parti- 
cle, is denoted ^ while the momentum of the associated (5*) hadron (typically a 
7r°) is denoted plj^ g. The azimuthal angle between the photon and the 77° is A^P2\ 



h^h dN''^^'''"Tl-'(J'i'i^ 



Pt.l 



(.55) 



The integrals are carried out within momentum intervals defined by p7p ^ and g. 
The yields in the above expression are defined in Eqs. ((5T|) and (|53l) . The correla- 
tion defined in Eq. ((54|) can be considered as a snapshot of the integrand in the 
coincidence probability defined in Eq. (j55p . In the same fashion, one can define a 
hadron-triggered coincidence probability.'^^ If the tt" is the trigger particle, the mo- 
menta are instead denoted by p7p g and which then become the lower limits on 
the integrals in Eq. (j55p . The away-side coincidence probability for the azimuthal 
correlation of photon-hadron pairs can also have a double-peak or single-peak struc- 
ture, depending on the trigger particle selection and kinematics.'^^ 
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Fig. 28. The 7 - 7r° correlation P(A0) defined in Eq. ^ in minimum bias p + p and p-l-Pb 
collisions at = 5 TeV and 'r]^ = rj-y = 3 (with the convention that the proton beam moves 
toward forward rapidity) obtained via the rcBK evolution equation with several initial saturation 
scales, Qq^ = ^Qop with TV = 3, 5, and 7, in two difi'erent transverse momentum regions: < 
(left) and p'^ > p^ (right). Taken from Ref.llHlEIZl 



In Fig. [221 predictions of the azimuthal correlations between the produced 
prompt photon and hadron, calculated employing the coincidence probability CPh 
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in p+Pb collisions are shown for = 0.2, 5, and 8.8 TeV. Equation (|55|) is used in 
the denominator of Eq. (j55p . The colhnear divergence was removed from Eq. (j55p by 
introducing the photon fragmentation function. The numerator is calculated using 
Eq. (j5ip . An overall normalization problem may result given that the away-side cor- 
relation at A(/) = 7r/2 is not sensitive to the colhnear singularity. Proper treatment 
of this problem requires a full NLO calculation which is currently unavailable. How- 
ever, choosing a different photon fragmentation functiorPslEiniElI] ^\]\ only slightly 
change the results given the freedom to choose the fragmentation scale and the 
rather large uncertainties due to Qq^. This possible normalization problem is not 
present in the correlation defined via Eq. (j54l) . 
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Fig. 29. The 7 — 77'' coincident probability, CP;, (Ai/i), defined in Eq. II55I I in minimum- bias p + A 
collisions at r?^ 3 (with the convention that the p roton beam moves toward forward 

rapidity) for = 0.2, 5 and 8.8 TeV. Taken from Ref.lMUll 

Although there are theoretical uncertainties associated with the strength of the 
photon-hadron correlations, due to both higher-order corrections and the less con- 
strained saturation scale Nevertheless, the decorrelation of away-side photon- 
hadron production increases with energy, rapidity and density. This decorrelation, 
together with the appearance of a double or single-peak structure, are robust pre- 
dictions of CGC/saturation effects in the leading-log approximation. 

6.4. Proton-nucleus dilepton and photon production at the LHC: 
gluon saturation and shadowing (R. Baier, F. Gelis, A. H. 
Mueller and D. Schiff) 

Here, a qualitative discussion on saturation and shadowing effects on photon and 
dilepton production at very large rapidity in high-energy proton-nucleus collisions 
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is given. The aim is to describe the various steps of assumptions leading to the 
present understanding. 

I'^ y/^NN — 200 GeV d+Au data on high-p^ hadron production at large rapid- 
ity (on the deuteron side) there is a significant suppression of hadron production 
compared to the expectation from p + p colhsions. This resuh suggests that there 
may be a significant amount of leading-twist gluon shadowing in the nuclear wave- 
function in the region probed by forward hard scattering at RHIC, with even more 
expected at the LHC. 

In many ways hard photons (or dileptons arising from virtual photons) are a bet- 
ter probe than high-p^ hadrons. Photons are less sensitive to fragmentation effects 
while final-state effects are almost absent. Thus, at ~ 2 — 3 GeV/c, leading- 
twist factorization is expected to be accurate. Hence the nuclear gluon distribution 
probes x values somewhat smaller than 10~^ in p-|-Pb collisions at the LHC. 

Increasing the photon rapidity into the forward (proton rapidity) region, ?/ > 0, 
the gluon x probed decreases rapidly, e.g. for virtual photons of mass M and rapidity 
y produced in the process g-f .g — 7> 7*, x ~ (M/y^) exp (— y). With ^/s = 5500 GeV, 
M — 5.0 GeV/c^, and y — 3.5, x — 2.7 x 10^^, so that indeed a fast quark in the 
direction of proton rapidity produces virtual photons that probe the small x gluon 
distribution in the nucleus. 

It is extremely interesting to explore, at least qualitatively, the size of the sup- 
pression one might expect in such reactions. The treatment here follows Ref.li^ 
which can be examined for further details. This discussion is based on a picture 
derived from the Color Glass Condensate (CGC) effective theory for the gluon dis- 
tributior0 together with BFKL evolution (see e.g. Ref.l^^ to reach higher values 
of y. 



6.4.1. Factorized formula for the inclusive 7* cros.s .section 

A dimensionless observable, a{Q,Y), can be written in the fcj^-factorizcd form 

aifi, Y) = Hiqr, ft) c^aiqT. Y) , (56) 

where fi is the hard scale of the reaction, equal to the transverse momentum, fcy, 
of the (massive) photon, i.e. 

.(..,y)..fe..,y).-pj-|j-^, (57) 

For simplicity incident quarks are considered instead of protons, q + A — )• 7* AT, see 
Fig. ISni Here kj- is the transverse momentum of the 7* , and z is the longitudinal 
momentum fraction of the 7* with respect to the incident quark momentum, z = 
/p+ , where p+ — )• 00 . The impact parameter of the q + A collisions is denoted b 
and Y = In 1/x. 



"'For recent reviews, see Refs.lHHlC^ll and references therein. 
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Fig. 30. Two gluon exchange graph for quark (p) + nucleus (P) production of a real or virtual 
photon, — 7(*)(fc) X. 

The leading-twist contribution involves only two exchanged gluons, as shown 
in Fig. 1301 H is the hard part of the interaction in the fc^-factorized form while 
(f'oiqT, Y) /q^ is considered to be proportional to the differential high energy q+A — > 
q + A cross section. 

6.4.2. Anomalous scaling and shadowing - specific predictions 

While there are theoretical uncertainties, there are also robust characteristic quali- 
tative results. The function (pGi^T, Y) can be approximated by the scaling function, 

0G(fcT, Y) = 0G(fcT/gs(6, Y)) « (kl/Ql{b, y)) ^""^ , (58) 

with anomalous dimension Aq — 0.37 and saturation scale Qs- In 6 = collisions at 
fixed kT and Y , an anomalous A dependence is predicted, together with shadow- 
ing/suppression of the ratio RpA, 

In order to obtain results at large photon rapidities based on the BFKL evolution 
in the presence of saturation, the rapidity dependence of the saturation scale is 
chosen to be compatible with phenomenology, 

Q2(6 = 0,y = y) = g^(r = 0)exp(Ay), A ~ 0.3 . (60) 

There is significant suppression expected, RpA — 0.5, especially when the 7^*-' is pro- 
duced at forward (in the proton direction) rapidities, e.g. y > 3 and fcy > 2 GeV/c. 
More detailed predictions of forward inclusive prompt photon production and semi- 
inclusive photon-hadron correlations in high energy proton-nucleus collisions at the 
LHC using the CGC formalism are presented in Refs.'^^'^^ 

Measurements of photons and dileptons, in addition to hadrons with moderate 
transverse momenta but at large rapidities in p-|-Pb collisions at the LHC will 
provide important information supporting the saturation picture of high density 
gluon dynamics at high energies and small values of x. 
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7. Jets 

In this section, several results for jet production and modification in media are 
presented. Multi-jet production is discussed and predictions for the jet yields are 
shown in Sec. 17.11 Cold nuclear matter effects on jet and dijet production are 
presented in Sec. 17.21 The angular decorrelation of dijets due to saturation effects 
is discussed in Sec. 17.31 Predictions of long-range near-side azimuthal collimation 
event with high multiplicity are given in Sec. 17.41 Finally, enhanced transverse 
momentum broadening in medium is described in Sec. 17.51 

7.1. Multi-jet Production (N. Armesto) 

The study of jet production in Pb+Pb coUisions at the LHCP^lHIifil and their 
medium modifications, are important for probing the properties of the hot and 
dense matter formed in heavy-ion collisions and is thus a very hot topic in heavy- 
ion physics. Therefore, studies of (multi-)jet production in p + A collisions are of 
great importance as a cold QCD matter benchmark. Here the jet rates in minimum 
bias p+Ph colhsions at the LHC (4-1-1.58 TeV per nucleon) are computed using the 
Monte Carlo code in Refs.li^^lli^ that implements fixed-order perturbative QCD up 
to next-lo-leading order. This code produces at most three jets and contains neither 
parton cascades nor hadronization corrections. 

In the computation, the renormalization and factorization scales have been 
set equal and fixed to fi = fip = M-R = ^t, where Et is the total trans- 
verse energy in the event. The central set of the NLO MSTW2008 parton den- 
sities (MSTW2008nlo68cl)ii^ have been used for the unmodified nucleons. The 
nuclear modification of parton densities was examined using the EKS9#2lll50] ^nd 
EPSOgNLC^Sl sets. The standard Hessian method was employed to estimate the 
uncertainties coming from nuclear parton densities. The precision of the computa- 
tion, limited by CPU time, gives statistical uncertainties smaller than 10% for the 
bin with the highest Et in the results shown here. The anti-fcT jet finding algo- 
rithm with R = O.d^^ was used. Only jets with Et > 20 GeV are considered. The 
uncertainties due to the choice of nucleon parton densities, isospin corrections, and 
scale fixing, together with the infiuence of the jet-finding algorithm and the choice 
of nuclear targets and collision energies, were discussed elsewher^i^li^ and are not 
considered here. 

Figure I3l1 shows the sum of the 1-, 2- and 3-jet yields, as well as the individual 
2- and 3-jet yields within four pseudorapidity windows (two central, one backward 
and one forward) in the lab frame as a function of the Et of the hardest jet within 
the acceptance, -Et hardest- The yields, computed for luminosity C — 0.25 x 10^® 
cm~^s^^ integrated over a one month (10^ s) run (corresponding to an integrated 
luminosity of 25 nb^^),^^'* are quite large. The yields are given on the right-hand 
vertical axes in Fig. 1311 High yields are required to study cold matter effects in high 
Pt multi-jet production. Figure l311 demonstrates that such studies are feasible. For 
example, in the backward region, —4.75 < rj < —3, sufficiently high rates for the 
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sum per GeV of H-2+3-jet events can be achieved for -Erhardest < 50 GeV. The 
only regions where the yields from certain channels are too low for statistically 
significant results are those for 3-jet events in the forward region and 2- and 3-jets 
in the backward region (only the yields for 2-jet events are shown). Note that the 
turnover in the 2-jet yields when i^Thardoat ~^ 20 GeV comes from singularities that 
appear in the NLO calculation due to large logarithms because Et^ ^ Et^ ) ■ 

Resummation techniques are required to obtain reliable results in this kinematic 
region. 

Nuclear modifications of the PDFs are very small, maximum 0(20%), hardly 
visible in the yields in Fig. |3T] Note that the wide uncertainty bands in the largest 
-^Thardeat ^Ihs lu thc forward and backward regions are due to statistical fluctuations 
in the Monte Carlo and do not have a physical origin. The corresponding hot nuclear 
matter effects in Pb+Pb collisions are expected to be much larger. 

Further work to compare these yields with experimental data would require 
consideration of hadronization corrections and the effects of background subtraction 
such as in the well-established jet area method-f^^ These should be the subject of 
future studies. 

7.2. Cold matter effects on jet and dijet production (Y. He, B.-W. 
Zhang and E. Wang) 

Inclusive NLCliSSHlSiJ jet and dijet production are studied in p-|-Pb collisions. Cold 
nuclear matter effects are included by employing the EPS09,^* DSll^^'^^*^'^ and 
HKNO'i^i^ parametrizations of nuclear parton densities. The numerical results for 
the inclusive jet spectrum, the dijet transverse energy and mass spectra, and the 
dijet triply-differential cross sections, all using a jet cone size R = 0.4 in minimum 
bias p-|-Pb collisions are shown. 

The left-hand side of Fig. |32] illustrates the inclusive jet spectra scaled by A^coii 
including the nuclear modifications in the central rapidity region, \y\ < 1. Results 
for i?ppb are also shown. The spectra with EPS09 and IIKN07 show an enhancement 
in the transverse energy range 30 < Et < 200 GeV. However, the CNM effects are 
negligible over the entire region of jet Et with DSll. 

The right-hand side of Fig. l32]displavs the rescaled dijet Et spectra as a function 
of Et2 with jet 1 at a fixed transverse energy of Eti = 100 GeV employing the 
same modifications of the parton densities as before. Both jets are within the central 
rapidity region, |yi| < 1 and \y2\ < 1. The cross sections in the range 10 < Et < 
200 GeV exhibit a peak near Et = 100 GeV. Again the spectra with EPS09 and 
IIKN07 are enhanced. On the other hand, the DSll set shows very small suppression 
and enhancement as shown in Fig 1321 

Figure |33] gives the rescaled dijet invariant mass Mjj spectra and the corre- 
sponding nuclear modification factors. Here Mjj is defined as [(X^-Pn)^]"'^^^ where 
the sum is over all particles in the two jets.'^^ The maximum rapidity of the two 
leading jets, lyjmax — inax(|?/i|, |j/2|) is defined so that lyjmax < 1- Jets with energies 
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Fig. 31. The 1+2+3- (black), 2- (red) and 3-jet (green) cross sections as a function of tfie Et of 
the hardest jet within the acceptance. Results in different pseudorapidity windows in the labo- 
ratory frame computed for minimum bias p+Pb collisions are shown. The dashed lines are the 
results without nuclear modifications of the PDFs. Results with EPS09NLCE1] (solid) and EKS98 
(dotted]p3E^ are also shown. The bands correspond to uncertainties computed using the Hessian 
method for EPSOqI^H The right vertical axes gives the scale for the corresponding yields with an 
integrated luminosity of 25 nb~^. See the text for further details. 

greater than 40 GeV in the rapidity range |y| < 2.8 are selected. The CNM effects 
with EPS09 and HKN07 also enhance the Mjj spectra over a wide Mjj region 
while the effect with DSll is modest. 

Figure also shows the rescaled dijet triply differential cross section^^^ and 
their nuclear modifications. The momentum fractions xi and X2 are defined as 

i£jct ztEjct 

where the sums run over all particles in the jets. The results in Fig. [33] are shown 
for the rapidity difference between the jets, y■^, ~ 0.5(|?/i — 2/2!) where yi and 2/2 are 
the rapidities of the two leading jets. The results as a function of x show obvious 
deviations between nuclear modification factors which implies that the dijet triply- 
differential cross sections could be a good observable to distinguish between different 
shadowing parameterizations.'^^ 
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Fig. 32. The inclusive jot (loft) and dijot (right) results in p+Pb collisions at y's = 5 TeV. The 
dijet results, presented as a function of i5y2 given for fixed energy Eq^i = 100 GeV for jet 1. 
The rapidity acceptance for all jets is \y\ < 1. The Ex spectra are shown in the upper part of the 
plots while the nuclear modification factors are shown in the lower panels of the figures. 
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Fig. 33. (Left) The dijet invariant mass spectra in p+Pb collisions. (Right) The dijet triply differ- 
ential cross sections in p+Pb collisions. The nuclear modification factors, -Rppb, are also shown. 
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7.3. Angular decorrelation of dijets as a signature of gluon 
saturation (K. Kutak and S. Sapeta) 

Here predictions for the emergence of saturatiorP^ effects on dijet production 
in p+Ph scattering at the LHC. The resuhs are based on a study performed in 
j^gfgli69][i70] r^Yie saturation scale characterizes formation of a dense system of par- 
tons. There is growing evidence that the phenomenon of gluon saturation indeed 
occursfi^'i^ This calculation employs a high-energy factorization formalism which 
accounts for both the high energy scale of the scattering and the hard momentum 
scale pt provided by produced hard system?^^ Dijets separated by large rapidity 
gaps^^"* are considered here. More specifically, the focus is on a case where one jet 
is measured in the central rapidity region of the detectors while the other is at large 
rapidity, see Fig.jSH Such final states probe the parton density in one of the hadrons 
at low momentum fraction x while the other is at relatively large x. Predictions of 
azimuthal decorrelations of dijets production in p+Ph relative to p + p collisions 
are presented. 

forward jet 
Pi 

Pi 

central jet 



Fig. 34. Jet production in the forward (assuming the proton moves toward positive rapidity) region 
in hadron-hadron collisions. 

Assuming, without loss of generality, that xi ~ 1 and X2 1 (if the proton 
moves in the direction of forward rapidity), the cross section takes the form 

C?CT _ \ - PT1PT2 

dyidy2dpTidpT2d.A(j) ^^8TT'^{xiX2sy 

X Mag^cdXlfa/Aixi, P^) 0g/B(a^2, fcr) . \ 7 (62) 

t ~r Ocd 

where 

= Pti + PT2 + 2pTlPT2 COS A(/> , (63) 

and A(f> is the azimuthal distance between the outgoing partons. The squared matrix 
element Magged includes 2 — > 2 process with one off-shell initial-state gluon and 
three on-shell partons a, c, and d. The following partonic subprocesses contribute to 
the production of the dijet system: qg — > qg; gg qq; and gg gg\^^ The off-shell 
gluon in Eq. ([62|) is described by the unintegrated gluon density (t)g/B{x2, k^), a 
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solution of a nonlinear evolution equatiorP^li^ that depends on the longitudinal 
momentum fraction and the transverse momentum of the off-shell gluon, kj-. 
On the side of the on-shell parton, probed at high proton xi, the collinear parton 
density /q/aI^^Ij /^^) is appropriate. 
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Fig. 35. Ratio of differential cross sections for central-forward dijet production at ^/s = 5 TeV 
as functions of azimuthal distance between the jets A</> for three different jet pj- cuts. On the 
left, the selected rapidities correspond to the CMS detector while, on the right, the rapidity range 
appropriate for ALICE is shown. The calculations have been made with the convention that the 
proton beam moves toward forward rapidity. 



Figure |35] shows the ratios of the differential cross sections for central-forward 
dijet production in p-|-Pb relative to p -I- p as a function of the azimuthal distance 
between the jets, Ai/i. In the region At/) ^ tt, the gluon density is probed at small 
kx, where it is strongly suppressed due to nonlinear effects. The ratio is a signal 
of saturation that is sensitive to the enhancement of saturation effects afforded by 
the larger A of the Pb nucleus. 

The left-hand side of Fig.[35lshows the results for various jet pr cuts: pTjct > 15, 
25, and 35 GeV/c while the jet rapidities are restricted to positive values within 
the coverage of the central and forward CMS detectors. On the right-hand side 
of Fig. [35l lower pT cuts, pr > 10 and 20 GeV/c, are used and the rapidity is 
restricted to positive values within the central and larger rapidity range of the 
ALICE detector. The results show that, when the two jets are back-to-back, the 
cross section ratio is significantly smaller than one. This is a consequence of stronger 
gluon saturation in the Pb nucleus than in a proton. Thus the unintegrated gluon 
distribution in the region of small and medium kx is more suppressed in Pb than in 
protons, as shown in Fig. 1351 In addition to the effects shown here, the region near 
A(j) ~ TT is sensitive to Sudakov effects which further suppress the cross section. 
Therefore, refinement along those lines could be envisaged in the future. However, 
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the ratio shown here is less sensitive to these corrections. It is worth emphasizing 
that the suppression due to saturation predicted in Fig. [35] is strong and extends 
over a large enough range for experimental observation, even if the region near 
A(f) — TT could be further refined. 

7.4. Predictions for long-range near-side azimuthal collimation in 
high multiplicity p+Pb collisions (K. Dusling and R. 
Venugopalan ) 

Rapidity-separated dihadron correlation measurements can provide valuable insight 
into the gluon dynamics of the nuclear wavefunction. For example, it was recently 
realizecP^ that a novel "near-side" azimuthal collimation in high multiplicity p + p 
eventJi^ is a consequence of nonlinear gluon interactions and was found to be 
in excellent agreement with computations in the Color Glass Condensate (CGC) 
effective field theory .C^lisni 

The situation in p+Ph collisions is more attractive. The degree of near-side 
collimation is highly sensitive to the fine structure of the unintegrated gluon dis- 
tribution. It can be shown that the collimated yield is enhanced by (Q^^/Qs) in 
asymmetric collisions. Here Q^^ is the saturation scale probed in the lead nucleus 
and can be estimated to scale with the number of participants, Q^'^^ ~ ^partQs'^; 
where is the saturation scale for the proton. Therefore, for high multiplicity 
events where many nucleons in the lead nucleus participate, a significant near-side 
collimation is observed. Similar to the p + p data, calculations within the CGC 
effective theory are able to explain the data.l^^ 

This last study combined the physics of saturation and BFKL dynamics in order 
to understand the near-side collimation and the recoiling away-side jet. Both are 
required for a full understanding of the azimuthal structure of dihadron correlations. 
The per-trigger yields shown in Fig. [32] are well described by this framework for the 
symmetric p^'^ p^^°'^ windows where data are currently available. Predictions 
for the asymmetric p!^'^ 7^ p^^°'^ windows are also shown. 

7.5. Enhancement of transverse momentum broadening (H. Xing, 
Z.-B. Kang, I. Vitev and E. Wang) 

Both initial-state and final-state multiple scattering lead to acoplanarity, or mo- 
mentum imbalance of the two leading final-state particles. To quantify this effect, 
the transverse momentum imbalance qx is defined as: 

qT^PTl+PT2, (64) 

with the average squared transverse momentum imbalance 

Here, da/dVS is the differential cross section with the appropriate phase space 
factor, determined separately for each process. For example, h\ p + A collisions. 
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Fig. 36. Correlated yield 1 / Nx^i^(f N / dA<t> after ZYAM (zero yield at minimum, used to remove 
the A<I>-independent pedestal) as a function of |A(/>| integrated over 2 < |At;| < 4 for the most 
central multiplicity bin > 110. The CMS datJiSEl have currently only been provided for 

the diagonal components, p'rP^ ^ p?^^°'^ , of the correlation matrix. The theory curves are the 
result of adding the BFKL contribution responsible for the away-side jet and the "dipole"-like 
Glasma contribution. The solid black curve is the result for Qg^ = 0.504 GeV^ on N^^^^ = 14 and 
the dashed green is for Ql = 0.336 GeV^ on AfP'^.j = 16. 
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dps = dyidy2dp^ for photon+jet production and dVS — dyidy2dpTidpT2 for 
photon+hadron production. 

The enhancement of the transverse momentum imbalance (or nuclear broaden- 
ing) in h + A [h — p, 7*) collisions relative to h+ p collisions can be quantified by 
the difference: 

^{Qt) = {qT)hA ~ {qT)hp ■ (66) 

The broadening A{q^) is a result of multiple quark and gluon scattering and is a di- 
rect probe of the properties of the nuclear medium. Both initial-state and final-state 
multiple parton interactions are taken into account to calculate the nuclear broad- 
ening A(g|,) for photon+jet (photon+hadron) and heavy quark (heavy meson) pair 
production in p -\- A collisions. 

The nuclear broadening A{q^) can be calculated in perturbative QCD. A specific 
method based on double parton scattering has been discussed in detail in Refs.'^'^^ 
The derivation discussed here closely follows Ref.'^ The leading contribution to 
the nuclear broadening comes from double scattering: either in the initial-state 
or the final-state. The contributions from these diagrams in the covariant gauge 
can be calculated to obtain the following expression for the nuclear broadening of 
photon+jet production in p + A collisions: 



A(4) 



(67) 



Ea,b fa/p{x')fb/A{x)H^t,^^^i^, t, U 

Here Tj[^^(x) — T^^y]^(x) (or T^p^{x)) are twist-4 quark- gluon (or gluon-gluon) cor- 
relation functions associated with initial-state multiple scattering, defined a^SOllissi 

X \ {pA\F+{y2 )^,{Q)^+%(y-)F+'^{y^)\pA) , (68) 

X ^ {pA\F^+{y2 )F-+{Q)F\{y-)F+'^{y^)\pA) . (69) 

The corresponding twist-4 correlation functions associated with final-state multiple 
scatter are t''JI{x) and 7j;f^(a;). They are defined as in Eqs. ((68|) and (|69)) . except 



the ^-functions are replaced such that,!^^^ 

6{y- - y-) 9{-y-) ^ 6{y- - y-) 6{y-) . (70) 

The broadening in photon + hadron production is calculated similarly. 

Nuclear broadening is evaluated in a formalism where multiple scattering con- 
tributes to the cross section via higher-twist matrix elements in the nuclear state. 
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This framework follows a well-established QCD factorization formalism for particle 
production in p + A collisions and has previously been used to describe cold nuclear 
matter effects such as energy loss, dynamical shadowing and broadening. This work 
differs from more generic parton broadening phenomenology in because the color 
and kinematic structures of the hard part are evaluated exactly. In particular, the 
nuclear enhancement of the transverse momentum imbalance is studied in dijet and 
photon+jet production. 




2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 
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Fig. 37. Nuclear broadening, A(q^), is shown for dijet (left) and photon+jet (right) production 
at specified rapidities in p + A collisions as a function of N^o\\ . Results are shown for j/i = j/2 = 2 
in = 5 TeV p+Pb collisions and j/i = 1/2 = 1 at y/s = 200 GeV d+Au collisions. (The 
calculations have been made with the convention that the proton beam moves toward forward 
rapidity.) The jet transverse momentum is integrated over 30 < px < 40 GeV/c at y/s = 5 TeV 
and 15 < Pt < 25 GeV/c at ^/s = 200 GeV. The red line shows the result for RHIC kinematics 
with scattering parameter = 0.12 GeV^ while the yellow band represents the variation of in 
the LHC kinematics. 

The results are shown in Fig.[37l The line gives the baseline ^{qj) determined 
from RHIC as a function of A^coii- The band is the result in the LHC kinematics 
for a plausible range of the scattering parameter , defined in Eq. (|29|) . The band 
is broader for dijet production and ^{(It) is also larger. There is also a somewhat 
larger deviation of the bottom edge of the band from the RHIC result for dijet 
production. 

8. Gauge Bosons 

Gauge boson production is discussed in this section. Sect ion [8 . II makes predictions 
for W and Z'^ pT and rapidity distributions in p + p and p+Ph collisions. The 
charge asymmetry as a function of the decay lepton rapidity is also shown. 
Section [8 . 2 1 describes a calculation of the resummed pT distribution while trans- 
verse broadening of vector boson production is discussed in Sec. 18.31 
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8.1. W and Z production and charge asymmetry (P. Ru, E. 

Wang, B.-W. Zhang and W.-N. Zhang) 

Production of the gauge bosons, ^W" and is discussed. The transverse 
momentum and the rapidity distributions in min-bias p+Pb colhsions at ^/s — 
5 TeV are calculated at next-to-leading order and next-to-next-to-leading order. 
The fiducial cross sections crfid are calculated in the fiducial phase space for vector 
boson production used by ATLAS Collaboration EHEm in both p-|-Pb and p + p 
collisions. 

The fiducial Z° cross section is the inclusive cross section p -\- Pb — >■ j^* -|- X 
multiplied by the branching ratio for Z^ /^* — > within the fiducial acceptance. 
Here X denotes the underlying event and the recoil system. The fiducial acceptance 
is assumed to be the same as that defined by ATLAS \n p + p coUisionJi^ with 
the following cuts on the lepton transverse momentum and pseudorapidity, and the 
dilepton invariant mass: plj, > 20 GeV/c; |?7'| < 2.4; and 66 < mu < 116 GeV/c^. 

The fiducial cross sections are the inclusive cross sections p+ Pb +X 
multiplied by the branching ratios for — > liy within the fiducial acceptance. Fol- 
lowing the ATLAS definition,'!^ the acceptance cuts on lepton and neutrino trans- 
verse momentum and pseudorapidity as well as the W transverse mass are: pip > 

20 GeV/c; |r/i,^| < 2.4; > 25 GcV/c; and mr = ^2pij,p^{l - cos{(f>^ - 0'')) > 
40 GeV/c2. 

Results are simulated to NLO (0(qs)) and NNLO (0(a^)) in the total cross sec- 
tion employing DYNNLD for Drell-Yan-like production in hadron-hadron collisions.!^ 
The MSTWNLO and MSTWNNLO proton parton densities are used, along with 
the EPSOtl^ and DSSZ^^^ shadowing parameterizations. 

Tables 18.11 and 18.11 shows the fiducial cross sections, normalized to their per 
nucleon values for better comparison. The NLO and NNLO values are shown for 
comparison, both for p-|-Pb and p + p collisions at the same energy. The ratio 
between the NNLO and NLO cross section, giving some indication of the theoretical 
uncertainty and the convergence of the perturbative expansion for gauge boson 
production, is ^ 1.02 showing that the higher order corrections are small. When the 
NLO and NNLO results are compared in p-|-Pb collisions with the EPS09 shadowing 
parameterization, a similar correction is found. The difference between the EPS09 
and DSSZ parameterizations, both calculated at NLO, is also quite small, on the 
order of 1%. There is a slight decrease in the total fiducial Z° cross section. There is 
a larger decrease for and an enhancement in the fiducial cross section in p+Ph 
collisions. This is less an effect of shadowing than it is of isospin since ud and 
ud — >■ and there are more d quarks in the lead nucleus, causing the enhanced 
cross section. 

Figure [38] shows the Z^ pT distributions in p+Vh collisions at both NLO and 
NNLO. The NNLO result is somewhat higher and not as smooth as the NLO 
calculations which appear independent of the choice of shadowing parameterization. 
Differences between the results with EPS09 and DSSZ, which can be attributed to 
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Table 5. The total vector boson production cross 
sections in the fiducial phase space, crfld, in units 
of nb in p + p collisions. The results at NLO and 
NNLO are compared. 





P + P o-fld (nb) 




Decay channel 


MSTWNNLO 


MSTWNLO 


Z e+e- 


0.339 


0.332 


W+ -s> e+u 


2.35 


2.30 


W~ — > e'l/ 


1.47 


1.44 



Table 6. The total vector boson production cross sections per 
nucleon in the fiducial phase space, cra^/ {Ni^i^) in p+Pb colli- 
sions. Columns 2 and 3 compare the results at NLO and NNLO 
calculated with EPS09 while columns 3 and 4 compare the 
EPS09 and DSSZ shadowing parameterizations at NLO. 





p+Pb afld/(Afcoll> (nb) 






MSTWNNLO 


MSTWNLO 


MSTWNLO 


Decay channel 


EPS09 


EPS09 


DSSZ 


Z^e+e- 


0.338 


0.328 


0.329 


W+ e+u 


2.06 


2.02 


2.05 


W~ — >• e~v 


1.54 


1.52 


1.53 



shadowing effects rather than isospin, are only apparent when the ratio RpPbiPT) 
is formed. The slight decrease in per nucleon yield in p+Pb relative to p -\- p seen 
in the total cross sections in Tables 18.11 and 18.11 are due to the lowest pT bin, 
Pt < 20 GeV/c. At higher pT, the ratio increases above unity. However, the effect 
is not significantly larger than 5% over the entire pT range. 

The individual and W~ pr distributions are shown at NLO and NNLO on 
the left-hand side of Fig. |39l Again, the NLO results are smoother than the NNLO 
calculations. The difference in the overall cross sections are clearly observable: in 
the fiducial range of the calculations, the cross section is w 35% greater than 
the W~ cross section. The ratio of p+Pb to p + p is shown on the right-hand side 
of Fig. |39l The W~ ratio is larger than unity and increasing strongly with py. The 
effect is due to the greater abundance of d quarks in the Pb nucleus (126 neutrons 
vs. 82 protons). The valence d quark distributions in the neutrons, equivalent to 
the valuence u quark distributions in the protons, have a larger density at relatively 
high a;, causing the observed increase. Conversely, the lower density of valence u 
quarks in the Pb nucleus causes RpPh < 1 over the entire px range. The sum of the 
two charged gauge bosons shows a trend very similar to that of the in Fig. [38l 
revealing a result closer to the true shadowing effect. Even though the W'^ cross 
section is greater than that of the W~ , the isospin effect on the W~ is large enough 
to make the ratio larger than unity at high pT- 

Figures l40l and |4T] show the and W~ rapidity distributions, normalized 

per nucleon. For comparison the p + p and p+Pb distributions are shown both 
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Fig. 38. (Left) Normalized differential cross section (l/crfl;j)(iio"fid/'^PT)- (Right) The suppres- 
sion factor Rpa{pt) in 20 GeV/c pT bins. 




Fig. 39. (Left) The normalized W = + W differential cross section (l/iTfid)(do"fld/(ipT)- 

(Right) The suppression factor RpAipr) for W production given in 20 GeV/c pT bins. 



at NLO and NNLO. Aside from numerical fluctuations at NNLO, the order of 
the calculation makes little difference in either the shape or the magnitude of the 
rapidity distributions. The p+p distributions are all symmetric around y ~ while 
the p+Pb distributions are peaked in the direction of the Pb nucleus. 

There are significant differences between the and W~ distributions even for 
p^-p collisions. The distribution is considerably broader with peaks away from 
2/ = 0, at \y\ ^ 1.7 due to the larger average momentum fraction x of the valence u 
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y 



Fig. 40. The Z*^ rapidity distribution (1/ {N^oii)){drjfi^/dy). Results are shown for both p+Pb and 
p + p collisions in the center-of-mass frame for both systems. 



quarks in the proton. The greater density of valence u quarks in the proton leads to 
the larger overall cross section. The isospin effect tends to make the W'^ rapidity 
distribution more symmetric since the valence u quark distribution in the neutron, 
equivalent to the valence d distribution in the proton, has a smaller average x and 
lower density which reduces the cross section while removing the peaks away from 
midrapidity. The W~ distribution in p + p collisions is both smaller and narrower 
than the W'^ . This distribution shows the strongest isospin effect with rapidity. 




Fig. 41. The W+ (left) and W (right) rapidity distributions {1/ {N coll)) {da f^^/dy). Results are 
shown for both p+Pb and p + p collisions in the center-of-mass frame for both systems. 
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These differences are reflected in the charge asymmetry, defined as {N-^+ — 
^w- ) I {Nw+ + ^w- ) s-nd shown as a function of the decay lepton pseudorapidity. 
The results are given in Fig. 1421 The p -\- p asymmetry is symmetric around y — Q 
with a strong dip at midrapidity. The origin is clear from the individual and W~ 
rapidity distributions. There is, however, a strong forward/backward asymmetry in 
p+Pb collisions. In the direction of the proton beam, the asymmetry follows that 
of the p + p result. It falls off and becomes negative in the direction of the lead 
beam. This trend is independent of the order of the calculation and the shadowing 
parameterization used. 



+ 
+ 

I 



1 ' I ' r 

MSTWNNLO 
MSTWNLO . 



p + Pb 




MSTWNNLO + EPS09 
MSTWNLO + EPS09 
MSTWNLO + DSSZ 



Fig. 42. The charge asymmetry, [Ny^r+ — Ny^r- )/[N^w+ ~^^w- )> a-s a function of the charged 
lepton pseudorapidity in both p + p and p+Pb coUisions in the center-of-mass frame for both 
systems. 



8.2. Nuclear modification of the transverse momentum spectrum 
of production (Z.-B. Kang and J.-W. Qiu) 

The Collins-Soper-Sterman formalisrrP^ is used to calculate production over 
the full pT range at the LHC, 

dcTA+B^zo+x 1 



dy dp^ (27r) 



dHe'P^-'' W{b, Mz, xi,X2) + Y(pT, Mz, xi, X2) .(71) 



The W term gives the dominant contribution when pr <C Mz while the Y term 
is perturbatively calculable, see Ref.,^®^ allowing a smooth transition from the re- 
summed low PT region to px Mz where the fixed-order perturbative QCD cal- 
culations work well. In Eq. (|7T|) . xi = Mzj \fs and X2 = e^^ Mzj \fs while W is 
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given b; 



W^{b,Mz,xi,X2) b<b 



^^^'^'^^''''''''^ ~ \ W''{bn..jMz\x,,X2)F^P{b,Mz,x,,x,;b^,^) b > b, ^^^^ 



-'max 



where 5niax ^1/ (few GeV) is a parameter that specifies the region in which is 
perturbatively valid, and F^^ is a nonperturbative function determining the large 
b behavior of W and is defined below. In Eq. ([7^ . W^{b,Mz,xi,X2) includes aU 
powers of large perturbative logarithms resummed from ln(l/fe^) to ln(M|)^®^ 

W''{b,Mz,xi,x2)=e-'^'^''^'^'>W''{b,c/b,xi,x2), (73) 
where c is a constant of order one,'^^'^^!! and 

fMi ,,,2 



Sib,Mz)= %■ 



M. 



\n[-A]A{a,(p))+B{a,{^i)) 



(74) 



2/f,2 

with perturbatively-calculated coefficients A{as) and B{as) given in RefP^ and 
references therein. The perturbative factor in Eq. (|73)) . {b,c/b,xi,X2), has no 
large logarithms. It is expressed as 

W^{b, c/b, xi, X2) ^ (70 ^ fi/Aixi,fi^c/b)fj/B{x2,fi = c/b) (75) 

i=q,q 

where co is the leading order qq — > Z'^ partonic cross section.'^Sll The functions fi/A 
and /i/B are the modified parton distributions given bjff^ 



ft/A{xi,n) / -pC,/a{xi/^,tJ-)(l)a/A{ttJ-) 

a "'^1 ^ 



(76) 



where is over a = q,q,g, (t>a/Ai^T fJ-) ^-re the normal proton or effective nuclear 
parton distribution functions (PDFs) and C'^/a = J2n=o^ijl('^s/T^)"' are perturba- 
tively calculable coefficient functions for finding a parton i from a parton a, given 
in Ref.ESSl 

The non-perturbative function F^^ in Eq. (17^ has the form. 



' b"^ 

F^P (5, Mz, a;i , X2 ; b^^) = exp << - In ' ^ "-'^ 



) [.91 {{h'r ~ ibLX) 



+ 52(6'-6Lx)]-52(&^-6Lx)|- (77) 

where the explicit logarithmic dependence, ln(M| &max/c^)' '^as derived by solving 
the CoUins-Soper equation. The (72 term is a result of adding a general power 
correction to the renormalization group equation while the g2 term represents the 
size of the intrinsic transverse momentum of active partons.^^^ 

The coefficients of the two terms proportional to 6^ in Eq. ((77|) can be com- 
bined,Q^ 

G2 = lnf— f^j 52+ff2, (78) 
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to sum the dynamical and intrinsic power corrections. By requiring the first and 
second derivatives of W to be continuous at 6 = f»max, the parameters a and gi 
in Eq. (j77p can be uniquely fixed, leaving only one parameter, G2, sensitive to the 
power corrections and other nonperturbative effects. Taking §2 = 0.25 ±0.05 GeV^ 
and 32 = 0.01 ± 0.005 GeV^, Gf = 0.324 GeV^. Predictions employing Eq. ^ 
are consistent with all p + p and p + p data from the Tevatron and the LHC.'^^ 

The EPS09 NLO parameterizatiorPSl is used to account for the leading-twist 
nuclear effects on the parton densities in p + ^ collisions. Following the method pro- 
posed in Rcf.,^^^ the nuclear-size-enhanced multiple scattering effects are accounted 
for by choosing g2 92 A^/^. Then for Z" production, Cf^ = 0.689 GeV^ 

In Fig.USl the predictions for Z° production are shown. The cross section includ- 
ing resummation in Eq. ()71|) is evaluated employing the CTEQ6M parton densities 
at factorization scale p, = Mt/2 = 0.5y^M^ + Pr- The Y term is calculated at 
NLO in as The upper panel shows the Z° production cross section as a func- 
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Fig. 43. boson production in p -|- p and p+Pb collisions at y'i = 5 TeV and y = 



tion of pt ■ The black dashed curve is the p + p baseline while the red solid curve 
shows the minimum bias p-|-Pb result. The blue dotted curve is the minimum bias 
p-|-Pb result without the A^^^ enhancement of 52 so that the nuclear-size-enhanced 
dynamical power corrections from multiple scattering are absent. The lower panel 
presents the nuclear modification factor RpA- 

In Fig. 1431 the red solid curves are almost indistinguishable from the dotted 
curves. Thus, power corrections are not important for Z'^ production at LHC ener- 
gies. Therefore, Z° production inp+A collisions is an ideal probe of the modification 
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of the parton densities in nuclei as well as of the high energy "isospin" effect. 

The cross section at y = is dominated by gluon-initiated subprocesses for 
Pt > 20 GeV/c. That is, RpA is an excellent observable to study nuclear modifi- 
cations of the gluon distribution, heretofore effectively unknown, especially at the 
values of x and /i probed by Z° production. At factorization scale /i = Mz, the 
EPS09 gluon shadowing factoil^ is less than unity (shadowing) for x < 0.005 and 
greater than unity (antishadowing) over a sufficiently large range, 0.005 < x < 0.2, 
as shown in Fig. |44l 





EPS09 














A=208 






— gluon 






— U-valence 









10"'* 10"^ 10'^ 10'' 1 

X 



Fig. 44. The EPS09 shadowing ratio at scale Q = Mz- Rf = ft/A(^,Q^)/fi/p(^,Q^)- The red 
solid curve shows the gluon ratio while the blue dashed curve shows the valence u-quark ratio. 



As indicated in the lower half of Fig. |43l the nuclear modification factor RpA 
for production is suppressed at pr < 10 GeV/c and enhanced at high pr, up 
to Pt ^ 70 GeV/c. The low pT suppression is a result of low x shadowing in 
EPS09. However, the strong enhancement over such a large pr range was unex- 
pected. After a careful examination of the kinematics, it was found that, at ?/ = 0, 
X ~ Mt/\/s ^ 0.013, already in the EPS09 antishadowing region. Thus, the clear 
enhancement of RpA in the large px region in Fig. 23] can be explained by the 
broad EPS09 antishadowing region at the Z^ scale. The large pr enhancement of 
RpA vanishes for other shadowing parameterizations with smaller gluon antishad- 
owing. Therefore, the measurement of RpA for production in p+Pb collisions at 
the LHC provides a clean and unique test of gluon antishadowing as proposed in 
the EPS09 parameterization. Furthermore, RpA is a direct measurement of nuclear 
gluon distribution since the cross section is dominated by gluon-initiated subpro- 
cesses for Pt > 20 GeV/c. 
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8.3. Transverse momentum broadening of vector boson production 
(Z.-B. Kang and J.-W. Qiu) 

Finally, transverse momentum broadening of inclusive vector boson production, 
A{pi) + B{p2) -> V[J/ip, T, Z°]{q) + X at the LHC is discussed. The average 
squared transverse momentum of vector boson production is 

(79) 



2 dOAB^V 



dy dq^ 

The transverse momentum broadening in p+Ph collisions is defined as 

^{qT)pPh{y) = {qT){y)pPh - {<it){v)pp ■ (80) 



Following the derivation in Refs-jliSSEMl the first nonvanishing contribution to 
the transverse momentum broadening of heavy quarkonium production is 

^/^2 \CEM ^ f ^T^'^as ^2 ^l/3\ (Cf + Ca) Oqq +2Ca <Jgg + AcTgg 

where the superscript "CEM" indicates that heavy quarkonium production is eval- 
uated in the Color Evaporation Model (CEM). A similar result was derived in 
the NRQCD approach.li^ The Uqq and Ugg partonic cross sections are contri- 
butions from quark-antiquark and gluon-gluon subprocesses, respectively.'^^ The 
Affgg term is a small, color-suppressed correction to the gg subprocess derived in 
j^gf llMI Jn the region where the gg subprocess dominates heavy quarkonium pro- 
duction, (Tgg ^ Uqq^ ^'^ggi heavy quarkonium broadening is further simplified ai^^ 

A(<z|)g™«2C^(^A^AV3^ . (82) 

In Fig. 1451 the predictions of transverse momentum broadening of Drell-Yan type 
vector boson production in p+Pb collisions at the LHC at y = are shown as a func- 
tion of iVcoii. To determine the effective dependence on A'coii, the A^l^ in Eq. (jSTj) is 
replaced by A^/^A'^coii(&)/^coii(&minbias)- In p-l-Pb collisions at the LHC, a Glauber- 
model calculation with crJJ\,^ — 70 mb at ^/s — 5 TeV gives coll {^min bias ) ~ 7. 
In addition to heavy quarkonium production, the broadening of W/Z^ production, 
calculated using the formalism derived in Ref . is also shown in Fig. |45l The 
dramatic difference in the magnitude of the broadening between heavy quarkonium 
and WjZ^ production in Fig. |45] should be a signature QCD prediction. 
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Fig. 45. The transverse momentum broadening of vector boson production in p+Pb collisions at 
y = 0, shown as a function of N^^ii . Th e T (red solid), J/tp (red dashed), (black solid), and 

(black dashed) results are given.EHS] 
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